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Abstract 

This article reviews developments in the theory of rapidly rotating degenerate 
atomic gases. The main focus is on the equilibrium properties of a single component 
atomic Bose gas, which (at least at rest) forms a Bose-Einstein condensate. Rotation 
leads to the formation of quantized vortices which order into a vortex array, in close 
analogy with the behaviour of superfluid helium. Under conditions of rapid rotation, 
when the vortex density becomes large, atomic Bose gases offer the possibility to ex- 
plore the physics of quantized vortices in novel parameter regimes. First, there is an 
interesting regime in which the vortices become sufficiently dense that their cores - as 
set by the healing length - start to overlap. In this regime, the theoretical description 
simplifies, allowing a reduction to single particle states in the lowest Landau level. Sec- 
ond, one can envisage entering a regime of very high vortex density, when the number 
of vortices becomes comparable to the number of particles in the gas. In this regime, 
theory predicts the appearance of a series of strongly correlated phases, which can be 
viewed as bosonic versions of fractional quantum Hall states. This article describes the 
equilibrium properties of rapidly rotating atomic Bose gases in both the mean-field and 
the strongly correlated regimes, and related theoretical developments for Bose gases in 
lattices, for multi-component Bose gases, and for atomic Fermi gases. The current 
experimental situation and outlook for the future are discussed in the light of these 
theoretical developments. 

Keywords: Bose-Einstein condensation; superfluidity; quantized vortices; fractional 
quantum Hall effect. 
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1 Introduction 

One of the most remarkable characteristics of a Bose-Einstein condensate (BEC) is its re- 
sponse to rotation. As was first understood in the context of superfluid helium-4[T], a 
Bose-Einstein condensate does not rotate in the manner of a conventional fluid, which un- 
dergoes rigid body rotation. Rather, the rotation of superfluid helium leads to the formation 
of an array of quantized vortex lines. Quantized vortices appear also in the superfluid states 
of helium-3 and in type-II superconductors in applied magnetic field0 which may be viewed 
as condensates of pairs of fermionic atoms or electrons. 

The experimental achievement of Bose-Einstein condensation in ultra-cold atomic gases - 
formed either from the condensation of atomic bosons [H [3], or of pairs of atomic fermions 4J 
- opens up a wide range of new features in the physics of quantized vortices and vortex 
arrays, allowing access to parameter regimes that are inaccessible in the helium superfluids 
or type-II superconductors. This can lead to novel properties of the rotating groundstates, 
including the possibility of exotic strongly correlated phases. In recent years, there have 
been advances in experimental capabilities and in the theoretical understanding of rotating 
ultra-cold atomic gases in these unconventional regimes. The aim of this article is to review 
the theoretical developments. It focuses mainly on situations in which there are large vortex 
arrays in a regime of rapid rotation, and on the theoretical predictions of the novel phases 
that can appear in this regime. 

We start with a brief introduction to the properties of Bose-Einstein condensates and to 
the physics of quantized vortices. 

1.1 Atomic Bose-Einstein Condensates 

Consider an ideal (non-interacting) gas in three dimensions, of identical particles of mass 
M, at number density n, and in equilibrium at a temperature T. The temperature sets the 
thermal de Broglie wavelength, At, via Ti 2 /MA T ~ fcsT, and the density the mean inter- 
particle spacing a ~ n -1 / 3 . At low temperatures, when At ~ a the gas must be described 
by quantum theory, and its properties depend strongly on the statistics of the particles. For 
bosons there is a phase transition when At ~ a, at a critical temperature 



For T < T c the gas is a Bose-Einstein condensate (BEC), characterized by a non-zero 
fraction of the particles occupying the same quantum state. 

1 The connection between a rotating neutral fluid and a charged fluid in a magnetic field will be clarified 
below. 




(1) 
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Until recently there was only one experimental realization of an atomic BEC. The tran- 
sition of helium-4 into a superfluid state below T c = 2.17K is known to be associated with 
Bose-Einstein condensation! 2 ] albeit in a system in which the inter-particle interactions are 
relatively large. However, in recent years, Bose-Einstein condensation has been achieved 
in a wide variety of atomic species. These are prepared as metastable low density gases, 
n ~ 10 12 — 10 15 cm -3 , confined in magnetic or optical traps. At such low densities, the BEC 
transition temperature is extremely small, T c ~ lOOnK. Nevertheless, by a combination of 
laser and evaporative cooling these low temperatures can be routinely achieved. 

At the low temperatures involved, the thermal de Broglie wavelength, At ~ a ~ 0.1 — 
1/im, is much larger than the typical range of the inter- atomic potential. The two-particle 
scattering is therefore dominated by s-wave scattering, with a scattering length a s that is 
typically of order a few nanometres (for 87 Rb, a s ~ 5nm). A typical atomic gas of bosons 
is weakly interacting, in the sense that na^ <C 1. Consequently, it can be well described as 
an ideal Bose gas with very small condensate depletion. (This is in contrast to superfluid 
helium, for which the strong interactions cause significant condensate depletion.) That said, 
inter-atomic interactions are still present and are important for many physical properties. 
Non-zero repulsive interactions will cause the gas to behave as a superfluid, with non-zero 
critical velocity [6]. This, in turn, leads to the appearance of quantized vortices when the 
superfluid is forced to rotate. 



1.2 Quantized Vortices 
1.2.1 Quantized Vortex Line 

It was noted by Onsager and Feynman that superfluid helium cannot rotate as a conventional 
fluid. A conventional fluid rotating at angular frequency ft has the velocity field of rigid 
body rotation 

v = n x r (2) 

for which the "vorticity" of the flow, V x v, is uniform 

Vxv = 2fl. (3) 

If, as is believed to be the case, the superfluid is described by a superfluid wavefunction 
ip s = y^nje"^*") then the superfluid velocity is 

« S = ^V0. (4) 

Hence, the fluid vorticity apparently vanishes: 

Vxd s = AvxV(^ = 0. (5) 

The last equality relies on cf) being a smooth function of position and overlooks the possibility 
that the phase <f> might have line-like singularities (point-like in 2D) around which cf> changes 
by an integer multiple of 2n. These are the quantized vortex HnesyQ. Integrating the vorticity 
over a 2D surface containing such a singularity gives 

V X v s ■ dS = <j> v s ■ dl = — x integer (6) 



2 Bose-Einstein condensation has been measured by neutron scattering, yielding a condensate fraction at 
low temperatures of about 9% [5]. 
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Figure 1: Schematic diagram of the velocity field around a quantized vortex line with one 
quantum of circulation, k = -h. The superfluid density falls to zero within the vortex core, 
which has a scale set by the healing length £. 



indicating a delta-function contribution to the fluid vorticity on the vortex line. The "cir- 
culation" of the vortex, defined as 

a ■ dl = -j^ x integer (7) 

is therefore quantized in units of -p. This leads to a characteristic velocity profile, with an 
azimuthal velocity \v 8 \ which diverges as r — > 0, see Fig. [TJ In order to avoid the infinite 
kinetic energy associated with this divergent velocity, in the core of the vortex the superfluid 
density n s tends to zero, over a lengthscale of the healing length, £. 

The healing length is an important characteristic of a superfluid. Within the Gross- 
Pitaevskii mean- field theory [6] for a homogeneous interacting Bose condensate the energy 
is 



E-fiN 



^IV^ + iglV-l 4 -^! 2 



d 3 r (8) 



where tp is the condensate wavefunction. This is written for an atomic Bose gas with s-wave 
interactions, with the contact interaction strength, g — Airh 2 a s /M, chosen to reproduce 
the s-wave scattering length a s [see Eqn. (fl"4)) ]. In many situations, Eqn.([8]) is an accurate 
microscopic description of an atomic BEC, owing to the small condensate depletion na 3 <C 1. 
For a uniform fluid, minimization with respect to \ip\ 2 gives fi — g\ip\ 2 . There is therefore 
one characteristic lengthscale of the equilibrium fluid, the healing length, set by 



h 2 , , , h 2 I 



2Me ^^ m = 7ma = s - (9) 

For an atomic gas, the healing length is typically of order £ ~ 0.5/im and is large compared 
to the typical particle separation, £/a ~ ^Jd/a s > 1. In the case of superfluid helium, for 
which interactions are strong, the healing length is very short, £ ~ 0.8A; the vortex cores 
therefore have a size of order the inter-particle spacing. 

1.2.2 Vortex Lattice 

If superfluid helium is forced to rotate, for example by cooling liquid helium in a rotating 
vessel from above the superfluid transition (where the conventional fluid rotates as a rigid 
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body) into the superfluid phase [7], the superfluid establishes an array of singly quantized 
vortex lines oriented along the rotation axis. 

For dilute vortices, spaced by a distance a v S> £, the arrangement of vortices is dominated 
by the influence of the kinetic energy of the superfluid flow. Minimizing the kinetic energy 
of the superfluid flow in the rotating frame leads to the conclusions [8] that: 

(i) The mean vortex density in the plane perpendicular to the rotation axis is set by 
the condition that the mean superfluid flow (averaged over a lengthscale large compared to 
the vortex spacing) mimics rigid body rotation. The integral of the vorticity of a uniformly 
rotating body ([2]) over an area A is 2 Aft; setting this equal to the circulation of a superfluid 
containing N v vortices in this area, N v h/m, leads to Feynman's result for the mean vortex 
density [T] 

_iV v 2Mil 
n v = — = — . (10) 

(ii) The vortices experience pairwise logarithmic repulsive interactions, with an interac- 
tion energy per unit length of 



£>(|Jii-iWO (ii) 



2ttM 

Kj 

where n s is the superfluid density and Ri = (JQ, YA is the position of the i th vortex in the 
plane perpendicular to the rotation axis. The groundstate configuration of a set of classical 
particles interacting with logarithmic repulsion is a triangular lattice. One therefore expects 
the vortices to arrange in a triangular lattice at low temperature [5]. 

The direct observation of vortex arrays in rotating superfluid helium poses a very sig- 
nificant challenge: it requires the measurement of a set of very small vortex cores, of size 
£ ~ 0.8A, spaced by the much larger lengthscale a v — 1/ y^V ~ 1mm (at the achievable rota- 
tion rates of superfluid helium) . Images of the positions of quantized vortices emerging from 
the surface of rotating helium have been obtained in very remarkable experiments [TTH HTj . 
showing arrays of small numbers of vortices. 



1.3 Rotating Atomic Gases 

Vortex lattices may be generated in atomic BECs by confining the gas in a magnetic trap 
of cylindrical symmetry, and stirring the gas with a rotating deformation |12[ 1 13L 1141 . The 
formation of vortices is driven by dynamical instability of the surface modes [15l [TBI Wf\ PT8] 
followed by a period of turbulent flow before steady state is reached. 

We shall focus on the equilibrium states that can be achieved, assuming that the system 
reaches a steady state at a rotation rate f2. Under typical conditions, the rotation rates are 
51 ~ 27r x 100Hz, setting a mean inter-vortex spacing of a v ~ 2/mt. This is small compared 
to the condensate size, allowing large numbers of vortices to be created in the condensate 
(observations of up to 160 vortices were reported in Ref.[13]). The gases are typically in 
the regime £ ^ a v for which the vortex lattice structure is determined by the kinetic energy 
of superfluid flow (|11|) and one expects a triangular lattice ordering at low temperatures. 
Indeed, images of the particle density, taken following expansion of the cloud, show large 
numbers of vortices ordered in a triangular lattice arrangement [T2l [T3l [T4l [T9l [20] . 

There are many reasons why cold atomic gases are very interesting systems in which 
to study the properties of superfluid vortices and vortex arrays. The aspects that we shall 
focus on in this review are: 

• High vortex density. In atomic BECs the inter-vortex spacing a v ~ 2/zm is compa- 
rable to the vortex core size, as set by the healing length which is typically £ ~ 0.5/im. 
One can therefore envisage |21j entering a novel regime in which the rotation frequency 
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is sufficiently high that the separation between the vortices a v becomes less than the 
healing length £. Indeed, this new regime has been achieved for a Rb condensate, by 
using an ingenious technique to achieve very high rotation rate and low density [20]. In 
superfluid helium, where the zero temperature healing length is £ ~ 0.8A, achieving 
this regime would require £1 > 10 13 rads _1 ! It is of interest to understand the properties 
of vortex lattices in the regime, a v <C £, which differ markedly from those of superfluid 
helium for which £ -C a v . 

• strongly correlated phases. One can envisage entering a regime of even higher vor- 
tex density, when the number of vortices is comparable to the number of particles 22 . 
In this case, theory predicts the formation of strongly correlated phases. These can be 
understood as bosonic analogues of the strongly correlated phases that are responsible 
for the fractional quantum Hall effect (FQHE) of electrons in semiconductors [23| [24] . 
These states have many interesting and novel properties. 

• Tunable interactions. There is the possibility to tune the inter-particle interactions: 
both the strength and the qualitative form of the interaction, by use of a Feshbach 
resonance and/or dipolar interactions. The groundstate of a rapidly rotating Bose gas 
is sensitive to the nature of the interaction. 

• Lattice potentials. The imposition of a rotating optical lattice 25, 26 , or the use 

of "artificial gauge fields" (27] EH [29] , allows studies of the interplay of vortex lattices 
with external lattice potentials. This raises interesting issues involving the effects of 
commensurability of the two lattice periods. Similar systems may also help to stabilize 
strongly correlated phases in experiment 29, 30, 31, 26 . 

• Multi-component (spinor) gases. The trapping and rotation of Bose condensates 
with multiple components is possible. This can give rise to more complex forms of 
vortex lattice order, including topological spin-textures, and raises the possibility of 
novel strongly correlated phases. 

• Fermi gases. Atomic Fermi gases can have interesting and unusual response to rota- 
tion. In particular, strongly interacting Fermi gases in the BEC-BCS regime behave 
as superfluids. Rotation leads to the formation of quantized vortices. Rapid rotation 
raises the possibility of interesting strongly correlated phases involving the interplay 
of the physics of superconductivity with the physics of the fractional quantum Hall 
effect (FQHE). 

As is clear from the above list of topics, this review will focus mainly on one aspect of 
the physics of rapidly rotating gases: the possibility of the appearance of strongly correlated 
equilibrium phases in novel parameter regimes accessible in experiments. Recent progress 
in this field has involved many other interesting developments. These include the non- 
equilibrium properties of rotating gases, dynamical excitations, the effects of confinement 
and trap geometry, and the many experimental developments. In the present article, these 
subjects will either be omitted or only touched upon briefly. The reader is referred to 
Rcfs.|4, 32J for excellent recent reviews that cover these topics in more detail. 

2 Rapidly Rotating Atomic Bose Gases 

The possibility of achieving high vortex density, and entering a regime where the mean 
vortex spacing a v is small compared to the healing length £ 21J, is one of the most striking 
differences between the physics of atomic BECs and that of rotating helium. We now describe 
the theoretical description of atomic gases in this rapid rotation regime. This discussion will 
form the basis for many of the theoretical results described in this review. 
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2.1 Rapid Rotation Limit 

The Hamiltonian describing an interacting atomic gas in an axisymmetric harmonic trap is 

H = E w + l M ^> + »?) + \ M "¥i + E V <r* (12) 



Most commonly, the interaction between ultra-cold atoms can be represented by a contact 
interaction]!)] 

V(r)=gS^(r) (13) 



with the interaction strength chosen as 

4nh 2 a R 



M 



(14) 



to reproduce the s-wave scattering length. However, in later sections we shall also consider 
situations involving dipolar interactions or a Feshbach resonance in the interaction. 

The Hamiltonian (|12[) describes N identical particles which are confined in a harmonic 
trap, with natural frequencies u± and u>\\ . Note that the potential is axisymmetric about 
the z-axis. We therefore (for now) neglect the influence of any deformation that is required 
to stir the gas. One can imagine that the system has been stirred, such that it has picked up 
a large angular momentum about the z-axis, but this stirring potential has been removed. 
Due to the axial symmetry, the angular momentum will be conserved. The question we 
shall address is: What is the groundstate of a Bose gas that has been prepared in this way? 
Equivalently: what is the lowest energy state of (|12p as a function of the angular momentum 
about the i-axis? 

For much of this review, we shall consider situations in which the interactions are rela- 
tively weak. In this limit, the relevant single particle states simplify. To motivate the form 
of the single particle states it is convenient to work not in terms of fixed angular momentum 
L, but in terms of its conjugate thermodynamic variable which is the rotation rate fl. 

2.1.1 Rotating Frame of Reference 

In a frame of reference rotating about the z-axis with angular momentum 1~2 = Qz, the 
Hamiltonian is [33] 

Hq = h — n ■ l . (15) 

In this frame, the one-body terms can be written in a suggestive way 

= ^+ 1 -M^+y^+ 1 -M^-n-rxp (16) 
\p-Mnxr\ 2 , 1 , 2 n2 . , 2 , a ,l 



, M/ + -M(u{-n 2 )(x 2 +y 2 ) + -Mu J 2 z 2 . (17) 

The kinetic term in this Hamiltonian is equivalent to that of a particle of charge q* experi- 
encing a magnetic field B* with 

q*B*=2MCt. (18) 

This connection shows that the Coriolis force in the rotating frame plays the same role as 
the Lorentz force on a charged particle in a uniform magnetic field 34] . Much use shall be 
made of the equivalence of these two effects in this review. Indeed, I shall frequently refer 
to "rotation" or "magnetic field" interchangeably, assuming that the connection is clear. 

One difference arises from the centrifugal force of the rotation, which has no analogue 
for a charged particle in a magnetic field. The centrifugal force acts to reduce the harmonic 
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confinement potential. Stability of a harmonically trapped gas requires that the confinement 
is retained in (fl7| , which requires that the rotation rate is below the "centrifugal limit" 

n<ujj_. (19) 

Analogies with the properties of a charged particle in a uniform magnetic field are most 
evident in the limit of a large number of vortices, N v ^> 1. In this case, very general 
considerations lead to the conclusion that f£ ~ ui± [see §2.5[ and Eqn. (|93|) ]. Essentially, 
one can then neglect any residual confinement and set f2 = uj±_ , viewing the system as being 
uniform in the plane perpendicular to the rotation axis. (The analogy holds also for finite N v 
with Q < luj_, but involves the Fock-Darwin spectrum for a charged particle in a magnetic 
field with harmonic confinement [35, 36 , which makes the present discussion somewhat less 
clear. This case is covered in the following section.) For il = oj±, the Hamiltonian (jTTJ) 
describes a quasi-2D system of particles in a uniform magnetic field. The energy spectrum 
takes the form familiar from studies of the QHE[2"3"| 

E n , m ,„n = 2tuo± ( n + i J + fiwii (nil + - J (20) 



2 J " V 11 2, 

where n = 0, 1, 2, . . . is the Landau level index, m — —n, —n + 1, — n + 2, . . . the angular 
momentum quantum number about the rotation axis, and nil = 0, 1, 2, . . . is the subband 
index for motion along the rotation axis. The spectrum is highly degenerate, with the single- 
particle states of different angular momentum m having the same energy. The effective 
"cyclotron frequency" is q*B*/M = 2uj±. 

The weak interaction limit occurs when the chemical potential /i ~ gh is small compared 
to the single-particle level spacing:^ 

fj, < , 2hu>±_ . (21) 

These conditions are equivalent to the conditions that the healing length is large compared 
to the inter-vortex spacing, £ 3> a v , and to the subband thickness in the z-direction, £ ^> an . 
Although (|21[) is typically not satisfied for a non-rotating gas, for a rapidly rotating gas the 
centrifugal forces spread the cloud out, the density falls and the system tends towards this 
weakly interacting regime [2D]. 

Under these conditions ([2Tj) Fl to a good approximation the single particle states are 
restricted to quasi-2D (ren = 0) and to the lowest Landau level (LLL, n = 0). The single 
particle wavefunctions are[2"T] 

i> m {r) CX (X + iy) m e -(*W)/2ai e -z 2 /2a\ ^ 

where 

s <23) 



3 For the strong inequalities in (1 2 1 I t to hold, the factor of 2 in this expression is clearly unimportant. We 
leave this factor in this formula merely as a reminder that the lowest energy single particle excitation out of 
the lowest Landau level has energy 2huij_. In practice, one expects the physics derived by theoretical studies 
in the regime ( 12 1 H to be at least qualitatively accurate even in experiments in which these inequalities are 
not well satisfied. 

4 The validity of the restriction of single particle states to the 2D LLL has been explored in Ref . 1371 . The 
numerical results are interpreted to indicate disagreement with the condition (12 1 1 for restriction to the LLL, 
for the case of incompressible states at fixed filling factor v = N/N v , Eqn J67l l. However, the criterion applied 
- the shift of the rotation frequency Q - appears to overlook the dependence of on Nv, e.g. Eqn, l|93 [l for 
./V v 3> 1. An analysis based on Eqn. H93H together with the assumption that fi ~ gn with a correction due to 
LL mixing 8fi ~ (gn) 2 /hu>± that is small (Sfi /i) leads to the conclusion that the quantity g ma x defined 
in Ref. 37 should vary as g max ~ N 2 /N. Thus g max oc 1/N at fixed AT V and g max <x N at fixed u, in rough 
agreement with the numerical results l37l . 
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are the trap lengths in the radial and axial directions. We refer to this regime (|21[) as the 
"2D LLL regime". 

To make connections to the quantum Hall effect as clear as possible, we shall introduce 
the complex representation 

x + iy 



c 

where I is the conventional magnetic length associated with q* B* — 2Muj± 

£ = 



2Mwi V2 ' 
The (normalized) 2D LLL basis states are then 



1 



V2vr2 m m! 



C m x 



-ICI74 



i (Trail) V* 



(24) 



(25) 



(26) 



For simplicity, in the following when writing wavefunctions in the 2D LLL, we shall omit the 
ubiquitous [bracketed] exponential terms, focusing only on the prefactors that are polynomial 
inC 



2.1.2 Laboratory Frame of Reference 

An alternative way to understand the restriction to the 2D LLL is to work in the laboratory 
frame, and consider the total angular momentum as the control parameter This ap- 
proach makes the appearance of Landau level wavefunctions somewhat less evident, but it 
has the advantage of clarifying that iV v ^ 1 is not required for the restriction to the states 
([22]) in the 2D LLL regime (|2"Tj) . 

For the harmonically trapped gas ()12|) the single particle energy spectrum is 

E = hujj_{2n± + \m\ + 1) + Uu\\{n\\ + 1/2) (27) 

where m = 0, ±1, ±2, ... is the angular momentum quantum number (about the z-axis) and 
n±, nil = 0, 1, 2, . . . are the radial and axial quantum numbers. For a set of particles i = 1, N 
occupying these states, the total angular momentum (in units of K) is 

N 

L = ^2rrn (28) 

i=l 

where m, is the angular momentum of particle i. 

For weak interactions (|21[) . to determine the groundstate at fixed total angular momen- 
tum L, we must first minimize the total kinetic and potential energies. Clearly, from [)27p. 
the particles then must only occupy single-particle states with n± = n\\ = 0. In addition, 
one can show that at fixed L (and assuming L > without loss of generality) particles must 
only occupy states with vm > 0[21]. The resulting set of single particle states - those states 
with n± = nil = and mi > - is precisely the set of 2D LLL states (|22|26|) . 

To see why one requires > 0, let us first write down the energy for a collection of 
N particles with total angular momentum ([28]) restricted to the 2D LLL (n± = nn = 0, 
Trii > 0). The total single-particle energy is 

N ( 1 \ 

E = ^2 [hu^im + 1) + ftu\\/2] = ( hLu ± + -Tiuj\\ J N + huj ± L . (29) 

i=l ^ ' 

Now, consider moving one particle, i = 1, from its state mi > to a new state m\ , which may 
have m[ < and then be outside of the 2D LLL. To conserve the total angular momentum 
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L, we shall, in addition, need to add mi — m'i units of angular momentum to the other 
particles. Keeping these particles in states with m > leads to an overall change in energy 

AE = huj ± [(-|mi| + |mi|) + (mi - mi)] = Huj ± [\m[ \ - m! x ] (30) 

which is an energy increase if m! Y < 0. Hence, at fixed L > the lowest energy states must 
have rrii > 0. 



2.1.3 Effects of Weak Interactions 

The inclusion of weak interactions poses a theoretical problem that is very closely related 
to that for electrons in the FQHE[23l [24]. For a system of N particles with a given total 
angular momentum, L, one must distribute the particles within the 2D LLL orbitals such 
that the angular momentum ([25]) is fixed. The single particle (kinetic+potential) energy 
of all these states (|29|) is the same, leading to a very high degeneracy at the single-particle 
level. The groundstate is determined by the action of the interactions within this degenerate 
space. This is a fundamentally non-perturbative problem, as in the FQHE. The differences 
are firstly in the nature of the inter-particle forces, and secondly, and most importantly, that 
here we are studying bosons, not fermions. We are interested in the nature of the phases 
that can emerge for rotating bosons for the specific forms of interactions that are physically 
relevant in ultra-cold atomic gases. 

For the most part, we shall focus on the case of contact interactions (fT3| between bosons. 
The interaction between particles is fully specified by the Haldane pseudo-potentials (see 
fJX|, which for contact interactions reduce to the single parameter 

This parameter sets the energy scale of the rotating bosons in the 2D LLL. 

For a set of N particles at total angular momentum L, we write the interaction energy 
of the groundstate of the contact interactions as Ej{L, N) = VqCi(L, N) where ei(L, N) is a 
dimcnsionless function. In the laboratory frame of reference, the total energy (single particle 
and interaction), is then, using Eqn.(|29|) 



E(L,N) = (hw_L + N + hu ± L + V ei{L,N) . (32) 

In the rotating frame of reference, the energy is 

Ea(L,N) = (hLu ± + ^hcj^ N + h{oj ± -n)L + V ei(L,N). (33) 

At given f2, the groundstate angular momentum L(Q) is obtained by minimizing (|33[) over 
the allowed integer values of L. Thus, in the rotating frame of reference, the groundstate 
angular momentum is a function of the dimensionless control parameter Tl(uj± — Q,)/Vq. From 
the viewpoint of the laboratory frame, in which the control parameter is angular momentum, 
the rotation frequency is a derived quantity, set by 

This is a standard thermodynamic relation stating that Q is the conjugate thermodynamic 
variable to L (expressed at zero temperature and in terms of the parameters we have defined) . 

For a highly anisotropic trap, with u>\\ <C there is the possibility to explore a related 
regime [38 

hujn < /j, < 2hu)j_ . (35) 
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The interactions are then small compared to the cyclotron splitting 2%uj±, so the motion per- 
pendicular to the rotation axis is still described by the LLL. However, since the interactions 
are large compared to the subband spacing, many subbands of the axial confinement are 
relevant. We shall refer to this regime (|35|) as the "3D LLL limit'' . Following the approach 
proposed in Ref . [38j . within mean- field theory in the 3D LLL regime, one can expect the 
density distribution along the rotation axis to adopt a Thomas-Fermi profile, with overall 
length 2W Z . As described in §2] the properties are then the same as in the 2D LLL, but 
with an overall change in interaction constant (|128|l . to V r TF = | M ^ a ^ ■ (See H2.6I for a 
discussion of the validity of mean-field theory in the 3D LLL regime.) 

In the above we have discussed the rapid rotation limit of atomic bosons. The same 
considerations apply equally well to rotating fermions, leading to 2D LLL and 3D LLL 
limits at sufficiently low density and weak interactions that (f2"Tj) or ([55)) apply. In §3] we 
shall describe the properties of rapidly rotating fermions in the LLL. 

2.2 Gross-Pitaevskii Mean-Field Theory 

The Gross-Pitaevskii mean-field theory for an interacting Bose gas amounts to the assump- 
tion that the many-body state is a pure condensate, in which all particles occupy the same 
single particle state ip(r) 

N 

*({ri}) = ]I^ r *)- (36) 

i=l 

As an approximation to the groundstate, the (normalized) condensate wavefunction ip is 
chosen to minimize the expectation value of the Hamiltonian. (We shall discuss the limits 
of applicability of this ansatz in §2.31 ) 

For a rapidly rotating Bose gas the single particle orbitals consist only of the states (|22p . 
Thus, the condensate wavefunction can be expanded in terms of these states 

ip(C,z) = ^2 c m ip m (37) 

m>0 

with \c m \ 2 = 1. The mean- field groundstate is obtained by choosing the coefficients c m 
to minimize the expectation value of the interaction energy, which for contact interactions 
(USD, is 

Ef p (L,N) = l -gN* J |^>(r)| 4 d 3 r (38) 
with a constraint on the (average) angular momentum, 

L = N^™\Crn\ 2 - (39) 

m 

The mean- field groundstate therefore depends only on L/N. The determination of the 
condensed LLL state that minimizes the interaction energy at fixed L/N constitutes the 
LLL mean-field theory. The application of GP mean-field theory within the 2D LLL regime 
(f2"Tj) , or the 3D LLL regime with TF profile (|3"5")) , is sometimes referred to as the "mean- 
field quantum Hall" regime in the literature. We avoid use of this terminology, preferring 
the term "mean-field LLL" regime. As discussed below, the physics of the mean-field LLL 
regime is in close analogy with that of the Abrikosov lattice in type-II superconductors 
(predating quantum Hall effects significantly!). We reserve the term "quantum Hall regime" 
to the regime of strong correlations, described in £12.41 where there are close analogies with 
the fractional quantum Hall effect. 

Note that, in general, Eqn. (|37[ ) is not an eigenstate of angular momentum, while the 
Hamiltonian conserves total angular momentum. The fact that the condensate wavefunction 
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Figure 2: Calculations of the particle density of an atomic BEC in the mean- field LLL 
regime, for two values of the angular momentum per particle L/N = 60 and 90, from 
Ref. [41] . The lattice structure is determined only by the interactions, which are here chosen 
to be contact interactions (Ti"5|) . Close to the centre of the trap, the vortices form a regular 
triangular lattice. [Reprinted figure from: N.R. Cooper, S. Komineas, and N. Read, Phys. 
Rev. A 70, 033604 (2004). Copyright (2004) by the American Physical Society. ] 

does not preserve the rotational symmetry of the Hamiltonian should not necessarily be 
viewed as a deficiency. Rather this wavefunction correctly captures the fact that the system 
spontaneously breaks rotational symmetry, in the limit N — > oo with L/N fixed. This will 
be discussed further in ^2.3.21 

Noting that the condensate wavefunction (|3T[) is a polynomial in £, it may be expressed 
in terms of its zeros [39 as 



where m max is introduced as a cut-off in the degree of the polynomial (|37p . and A is a 
normalization factor. These zeros are the complex positions of the quantized vortices. Thus 
the LLL Gross-Pitaevskii wavefunction is fully described by the positions of the vortices. The 
process of choosing the m max + 1 complex coefficients c m (with normalization) is equivalent 
to the choice of the locations of m max vortices. 

The minimization for contact interactions has been implemented numerically by several 
authors [40l |4"T1 [42] . A simple result is for the state at L/N = 1, for which the lowest energy 
condensate wavefunction is found to be 



in which all particles are condensed in the m = 1 orbital. For other values of L/N the 
system spontaneously breaks rotational symmetry. For large L/N the number of vortices 
grows, as N v ~ 3L/N, and forms a triangular vortex lattice that is weakly distorted by the 
confinement, see Fig. [2] Although small, the distortions of the triangular lattice are crucial 
to obtain the correct density distribution in the trap. Fixing the positions of the vortices 
on a triangular lattice leads to a density distribution that, when averaged over a lengthscale 
larger than the vortex spacing, is Gaussian[3S]. The true averaged density distribution is 
Thomas- Fermi like (i.e. an inverted parabola) |43j as is found from numerical studies [4Tj. 
The density distribution in this regime is discussed further in !j2.51 



1KC,Z)=A ]][ (C-Ca) 



(40) 



^{Q,z) = AQ 



(41) 
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The appearance of this triangular vortex lattice is very closely related to the appearance 
of a triangular lattice of flux lines in type-II superconductors close to the upper critical field. 
In that case, the order parameter is determined by minimizing the Ginzburg-Landau energy 
functional (per unit length L along the field direction) 



F-F n 



d 2 r . (42) 



Close to H C 2 the order parameter is small, so to a first approximation one can neglect 
the quartic term and minimize the quadratic terms. These describe electron pairs (hence 
charge — 2e) in a uniform magnetic field: the lowest-energy states are the (degenerate) 
lowest Landau level orbitals. As explained by Abrikosov[44] the groundstate is selected as 
the linear combination of the lowest Landau level orbitals that minimizes the quartic term. 
The solution is that the groundstate is a triangular vortex lattice [43] . and is characterized 
by the parameter 

^«H- 1596 ' <43) 

where (. . .) av denotes the (spatial) average over the unit cell of the lattice. This is mathe- 
matically equivalent to the LLL mean- field theory for rotating bosons, at least in the uniform 
case (i.e. close to the centre of a trap containing many vortices). The interaction energy 
([38]) is 

E? p = ±g 2d N 2 J\^(x,y)\Wr, (44) 

where g2d = 9 / (V^a\\) is the 2D coupling constant, and ?p(x,y) is 2D component of the 
condensate wavefunction, which is restricted to states in the LLL. Using the above result 
(|4"5|) , one concludes that the groundstate is a triangular lattice, and its energy (for a system 
of N particles in a total area A) is 

1 N 2 N 2 

Ef p = -g 2df 3 A — = (3 A V — . (45) 

In the last equality, we have replaced the area A by the number of vortices N v = Aj (ira 2 ^). 



2.2.1 Landau Level Mixing 

The limit of weak interactions (f2"Tj) where the LLL regime is applicable, corresponds to the 
case of dense vortices, £ ^> a±. For dilute vortices, ^<oi, interactions lead to contributions 
from higher Landau levels to the condensate wavefunction. We have explained above that 
for both (Cai (dilute vortices) and £ » a± (dense vortices) the mean-field groundstate 
is a triangular vortex lattice. As £,/a± varies between these two limits, there is a smooth 
crossover [4^1 147] . This crossover has been observed in experiments of Eric Cornell's group, 
which have achieved a chemical potential /i ~ gn that is less than the effective cyclotron 
splitting 2Hlo. In these experiments, /i remains of order 1.57iu)ii, so it is approaching the 
quasi 2D limit. Experimental evidence for the crossover into the LLL regime has been found 
in the apparent core size of the vortices 48], which shows the variation expected within 
calculations in the LLL [IS] H7J . The frequencies of the collective modes were predicted to 
show a reduction on entry into the LLL regime jlHl HO] ■ The experimental results appear 
to confirm the expected variation 51 , but this conclusion has been queried in a subsequent 
theoretical study that reassesses the expected shear modulus [52] and indicates that in the 
experiments of Ref. [5T| the collective mode frequencies have not reached the limiting values 
expected in the LLL regime. 

Very different physics controls the energetics in the two limits: for strong interactions 
£ < aj_ the vortex lattice is determined by the kinetic energy, for weak interactions in the 
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Figure 3: Geometry of two atoms with aligned dipole moments, as discussed in the text. 



LLL regime £ £ a± the vortex lattice is determined entirely by the interactions. One should 
therefore view the smooth crossover in the properties as somewhat of a coincidence. As 
discussed in ^2 .2.31 changing the form of the interaction potential can lead to large changes 
in the vortex lattice structure in the LLL regime, £ > aj_. 



2.2.2 Anisotropic Traps 

In order to model the process of stirring the trapped atomic BEC, one should include a 
potential that breaks the rotational symmetry of the Hamiltonian (fT2f . A simple form is 
to introduce an elliptical perturbation of the transverse confinement, which is at rest in the 
rotating frame of reference. This gives two (slightly) different oscillator frequencies u> x ,u> y 
in place of the single 

For this model, there exists a generalization of the LLL description for a weakly inter- 
acting gas 53J. The centrifugal limit (|19[) is now set by the condition that f2 < xnin.(u> x , u> y ). 
Close to this limit, the gas elongates in one of the directions transverse to the rotation axis 
leading to a quasi-lD geometry 54, 55, 56] [57] for rapid rotation. 



2.2.3 Dipolar Interactions 

Under certain circumstances, it has been proposed that the atoms (or molecules) in a cold 
trapped gas may carry an electric or magnetic dipole moment [58 . One should then add 
to the usual contact interaction the long-range dipole-dipole interaction. For two dipolar 
atoms with aligned dipole moments the effective interaction is 

V{r) = ^P(T)+C d ±=*%*° (46) 
M r a 

where 9 is the angle between the dipole moment and the line separating the atoms, see 
Fig.H 

A BEC of atoms with significant dipolar interactions has been realized by condensing 
chromium-52[59] , which is an atom with a very large magnetic dipole moment, /i = 6/ie, 

2 

such that Cd = % ■ The relative size of dipolar and contact interactions is parameterized 
by the dimensionless ratio 

e dd = —a— (47) 
3n a s 

which is edd — 0.16 for native conditions in 52 Cr, and has been further increased to e^d — 1 
by using a Feshbach resonance to reduce q c [60] . Very recently, there has been the exciting 
achievement of the condensation of hetero-nuclear molecules in their ground ro- vibrational 
state|61). Such systems can have large electric dipole moments, leading to very strong 
electric dipole interactions |58j. 
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Figure 4: Mean-field groundstates of a rotating BEC in the 2D LLL, for particles interacting 
by both contact and dipolar interactions. The relative size of these is controlled by a = 
V2/V0. The structure of the vortex lattice varies from (a) triangular (0 < a < 0.20) to 
(b) square (0.20 < a < 0.24) to (c) "stripe crystal", with simple rectangular unit cell, 
(0.24 < a < 0.60) to "bubble crystal" phases (a > 0.60) the simplest of which is shown in 
(d). [Reprinted figure from: N.R. Cooper, E.H. Rezayi, and S.H. Simon, Phys. Rev. Lett 
95, 200402 (2005). Copyright (2005) by the American Physical Society. ] 



Consider a rapidly rotating atomic gas (in the 2D LLL limit), and for simplicity let us 
choose the dipole moments to be directed parallel to the rotation axis [62]. The mean-field 
theory requires one to minimize the expectation value of the interaction energy 



\if>(r)\ 2 V(r - r')\^(r')\ 2 d 3 r d 3 r' (48) 



for ijj(r) m the lowest Landau level and for the interaction potential ([46]) . Although this is 
a simple generalization of the Abrikosov problem, it is one that does not naturally appear 
in that context, where the microscopic physics determining the superconductivity acts on 
scales much less than the vortex lattice period. 

The mean-field groundstates have been found for the case of a translationally invariant 
vortex lattice 62, 63J. The results[62j are shown in Fig. 2] as a function of 



where 1 1 



Vo 



h 2 a s / 2 Cd Itt Cd 



(49) 



Vo = V^]^ + V^^"V24 (50) 

piF (2m - 3)!! C d 

Vm>0 - V2 m!2™ 4 (51) 

are the Haldane pseudo-potentials for the dipolar interaction (|46[) in the limit a\\/a± <C 1. 
(See £jA]for definitions of the Haldane pseudo-potentials.) 

The results show that the mean-field groundstate is very sensitive to long-range repulsion, 
passing through a series of vortex lattice phases as a increases. The contact interaction 
aims to make |'0(^)| 2 as uniform as possible, while the long-range repulsion causes ^(r)! 2 
to cluster, leading to crystalline phases of clusters of particles at large a. These are referred 
to as "bubble crystal" phases [52], in analogy with the terminology used for structures of 
similar symmetry in 2D electron gases in high Landau levels pM] 165]. There is a sequence of 
bubble crystal phases as a increases, which are parameterized by the number of vortices per 
bubble, q. For 2D particle density ri2d the number of particles per bubble is B = n.2d-^a 2 , 
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B particles per "bubble" 



Figure 5: Schematic diagram of the bubble crystal phase. 



where a is the lattice constant of the bubble crystal, assumed triangular (see Fig. [5]) . For 
B ^> 1, the energy per particle is 



E 
N 



1. 



-v b 



-n 2 d 



V Q a 2 



KC d 



(52) 



where K — 11.034 is the Madelung constant for dipolar interactions on the triangular 

/ \ 1//3 

lattice[S5]. The energy is minimized for a m i n = ( ) , such that the number of vortices 
per bubble is 



1 V3 f KC A 



ira? ± 2 V 2Vb 



2/3 



= V3 



2/3 



(53) 



Vortex lattices have not, as yet, been created in experiments on the dipolar condensates. 
If these condensates can be made to enter the regime of rapid rotation that has been achieved 
in rubidium |51j one expects not just to see a crossover in the physical properties, but a true 
phase transition into a vortex lattice of different symmetry that depends on the ratio of 
short- to long-range interactions (j49| . The dependence of the mean- field phases on the 
degree of Landau level mixing, as quantified by /i/huj±, has been studied in Ref. |67j . 



2.3 Beyond Mean-Field Theory 

I now turn to describe the properties of rapidly rotating atomic BECs beyond the Gross- 
Pitaevskii mean-field approach. Much of our understanding has been derived from numerical 
exact diagonalization studies. These numerical techniques are described in Appendix [Bl 

In the remainder of this section I describe some of the results for rotating bosons in a 
harmonic trap (|12p in the weak interaction limit (|2ip , and concentrate on the case of contact 
repulsive interactions (|13|). 



2.3.1 Low Angular Momentum L < N 

For small values of the total angular momentum L the exact groundstates of a rapidly 
rotating Bose gas are known analytically. The following states (|54I55I57|) have been shown 
to be exact eigenstates of the contact repulsion by both analytical[68j and numerical [69] 
studies, and to be the groundstates for a class of repulsive interactions [70j [71] that includes 
contact repulsion (p~5|) . 
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• For L = there is only one state in the 2D LLL: the pure condensate in the m = 
orbital 



#L=o({Ci» oc 1. (54) 

• For L = 1 the groundstate is the centre-of-mass excitation of the L = state 

*L=i({Ci}) « Cc (55) 

Cc = (56) 

i=l 

for which the interaction energy is the same as for the L = groundstate, e/(l, N) = 
ei(0,N). 

• For 2 < L < N the groundstates are the elementary symmetric polynomials in the 
variables (d — Cc) 

« £ (Cfi-Cc)(c P2 -Cc)...(c M -Cc). (57) 

Pl<P2<Pi 

For contact interactions, the energies of the states with L < N (L ^ 1) are 68, 69\ [72] 
E X (L, N) = V ei(L, N) = ^V N(N - 1 - L/2) , (58) 

A close connection of these exact states for finite N to the condensed wavefunctions of 
the Gross-Pitaevskii theory can be made by considering the limit N — > oo. Consider first 
the case L = N, for which 

*i =JV ({Ci}) oc IJ(Ci - Cc) ■ (59) 

i 

If the co-ordinate C c were simply a number this would be the wavefunction of a pure conden- 
sate with a vortex at complex position £ c , Ean. (|40p . The fact that Cc is the centre-of-mass 
co-ordinate (|56[) . and therefore a function of the particle co-ordinates, means that this state 
is not fully condensed. Nevertheless, the fluctuations of the centre-of-mass from its average 
value (( c ) — 0, computed for the state ^l=n, are 

(C C 2 )~W. (60) 

Thus, in the limit N — > oo 

*l=n({Q}) N ^Y[Q (61) 

i 

which is the fully condensed state with a single vortex at the origin, (141 [) . This observation 
- that one recovers a fully condensed state in the limit N — > oo, keeping L/N = 1 fixed 
and staying within the 2D LLL (gfi <C ?ia;_i_, fi-wy) - is a specific example of the result that 
Gross-Pitaevskii theory is exact [73 in the limit N ~~ * oo with gn and L/N finite. One 
may generalize this to cases with L < N. To this end, we note that in the large N limit of 
(|54I55I57[) . obtained by setting ( c — * 0, the wavefunctions (|54I55I57[) are the set of all N + 1 
symmetric polynomials of {Q}. These polynomials can be generated by the condensate 
wavefunction 

*z({{i}) = l[(Ci-Z) (62) 

i 

where Z is a complex co-ordinate which may be viewed as the location of a single vortex. 
Thus, we conclude that, in the 2D LLL, the lowest energy condensate wavefunction for 
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L/N < 1 is precisely described as a single vortex state. For \Z\ finite and non-zero, the 
condensed state contains a single vortex located away from the origin, so it breaks rotational 
symmetry. 

The relation between the exact spectra and the condensed states has been studied in 
detail for L = 2N in Ref.[74]. 



2.3.2 Spontaneous Symmetry Breaking 

It is a general feature of the Gross-Pitaevskii theory for a rotating gas in an axisymmetric 
trap that the lowest-energy condensed state spontaneously breaks the rotation symmetry 
of the Hamiltonian. For example, for (|12| the general mean-field state is not rotationally 
invariant about the z axis[40], but shows small clusters of vortices with local crystalline 
order. 

Rotational symmetry breaking may be seen in the properties of the exact energy spectrum 
in the limit N — > oo with N v ~ L/N constant. Specifically, in this limit, the exact energy 
per particle (|33|) of the groundstate tends to the form 



En ^ N ^ const. + V Nf(L/N) + h(u ± - Q)(L/N) (63) 

for a certain subset of angular momentum states L, discussed further below. In this equation, 
Vo is the energy scale (|3"Tj) and f(L/N) = lirujv— ^e^i, N)/N 2 is a dimensionless function 
which is set by the form of the interactions. In the thermodynamic limit, we require that 
the energy per particle (and hence the chemical potential) remains finite. For example, to 
remain within the 2D LLL (|2ip requires the chemical potential to remain finite, /i < huj±, 
Tuvu. Since fi ~ VqN/N v , this requirement sets the condition that VqN — > constant. Under 
this condition, the function / has an overall scale of order 1. 

At a general value of CI, the exact groundstate angular momentum L* is found by min- 
imizing the energy per particle l|63p over the allowed integer values of L. Provided that 
f(L/N) is a (locally) smooth function such that the total energy has a simple minimum 
at L* - which is found to be the generic case - all states which have L that is within of 
order y/N of the optimal value L* have an excitation energy per particle that is of order 
Vo, and so vanishes in the appropriate thermodynamic limit, N — > oo with VqN — > const.. 
That is, in the thermodynamic limit there is a degeneracy between of order \/N states of 
angular momentum close the optimal value L* . The emergence of a degeneracy between a 
large number of states of different angular momenta is a signature of spontaneous rotational 
symmetry breaking in the thermodynamic limit [75] . 

Evidence for the appearance of rotational symmetry breaking can be found in the exact 
spectra at L > N, computed numerically for small numbers of particles. 

Fig.[6]shows the exact groundstate energy as a function of total angular momentum L for 
N = 25 particles, at a value of CI chosen such that the mean-field state is a pair of vortices, 
with a two-fold rotational symmetry. A notable feature is that the groundstate degeneracies 
occur between states which differ in angular momentum by multiples of AL = 2. This 
feature is an indication that in the thermodynamic limit the groundstate does not break 
rotational symmetry completely, but retains a two-fold symmetry. This is consistent with 
the mean-field state (see Ref.[10] at L/N — 1.75). For a state that retains a two-fold rotation 
symmetry, the values of L for which the limiting expression f|63|) applies are those that are 
spaced by AL = 2 from the true groundstate value L*; the function f(L/N) describes the 
smooth envelope of the dependence of the energies of the states at these angular momenta. 

The appearance of the quasi-degeneracies with increasing N is conveniently analysed 
by constructing scaling plots of [E(L + 1) + E(L — 1) — 2E(L)], as a function of the scaled 
angular momentum, L/N, see Fig.[7J This plot shows regions with quasi-degeneracies spaced 
by AL = 2 (1.5 < L/N < 2) and by AL = 3 (2 < L/N < 2.5), which indicate rotational 
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Figure 6: Groundstate energy as a function of the angular momentum, computed from exact 
diagonalization studies of rotating bosons in the 2D LLL with contact repulsion (|13|). at 
rotation frequency Q set to a value for which the mean-field groundstate shows two vortices. 
The groundstate is at L* = 44, the spectrum shows the appearance of low-lying excited states 
at angular momenta that differ from this by multiples of 2. In the thermodynamic limit, 
N — > oo, the exact spectrum shows the appearance of quasi-degeneracies between states 
spaced in angular momentum by AL = 2. This is a signature that in the thermodynamic 
limit the groundstate breaks rotational symmetry partially, leaving a two-fold symmetry 
axis [75], 

symmetry breaking to states with residual 2-fold and 3-fold rotation axes. Analysis of the 
second derivatives [E(L + AL) — E(L — AL) — 2E(L)] for AL = 2,3 in these regions (see 
Fig. [7] for the case AL = 2) allows the construction of scaling plots that confirm the above 
behaviour of the function f(L/N) at large N. 

For L/N < 1, rotational symmetry breaking can be deduced directly from considerations 
of the exact analytic wavefunctions of £12.3.11 For L < N, the exact energy spectrum for 
contact interactions ([55)) leads to f(x) = —x/i. Since the function f(x) has vanishing second 
derivative, this case is special from the point of view of rotational symmetry breaking. At one 
special value of the rotation frequency, 1 — fl/u>± = V$N / (4huj±), all states with L < N are 
degenerate in the thermodynamic limit. There are therefore ~ N degenerate states, rather 
than ~ VN. This enhanced degeneracy reflects the fact that, at this rotation frequency, the 
mean field state (|52"j) is itself degenerate with respect to the position of the vortex, that is 
with respect to all possible angular momenta < L/N < 1. 

Evidence of rotational symmetry breaking can also be found in the exact spectrum of 
the excitations above the lowest energy state at a given L. The relationship of the exact 
spectrum to the expected properties of Bogoliubov excitations of the mean-field groundstate 
has been discussed for L < N in Ref . [77] . 

The condensate fractions and condensed wavefunctions for small numbers of rotating 
bosons have been investigated in Ref .[78 . These are deduced from the exact groundstates, 
allowing for rotational symmetry breaking. The positions of the vortices in the condensate 
wavefunction (|40|) thus obtained show patterns similar to those from mean-field theory [40], 
The condensate fraction is large for a single vortex (L — N) for N > 9. For the system 
size studied (N = 6) the condensate fraction falls quickly to small values as the angular 
momentum is increased beyond the case of a single vortex (L > N) [78] . 
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Figure 7: Scaling plots of the (numerical second derivatives of) the groundstate energy 
as a function of the angular momentum, computed from exact diagonalization studies of 
rotating bosons in the 2D LLL with contact repulsion (fT3|) . The data collapse indicates 
that the large N asymptotic behaviour has been reached. These plots determine the scaling 
of the energies with N discussed in the text, which lead to the prediction of spontaneous 
rotational symmetry breaking in the thermodynamic limit N — > oo, VqN — > constant. There 
are regions which show clear two-fold and three-fold periods, corresponding to mean-field 
groundstates that break rotational symmetry partially to retain two-fold and three-fold 
symmetry axes respectively [76]. 



21 



Anisotropic Traps 

For a rotating trap with a quadrupolar deformation 54, 55, 56, 157[ 153] . which can rep- 
resent the stirred atomic BEC, angular momentum is not conserved. However, the wave- 
function parity is preserved: the deformation does not couple states with odd and even L. 
Close to the rotation frequency at which the first vortex enters the system, the exact spectra 
show close degeneracy of levels of opposite parity for large numbers of particles [79]. This 
quasi-degeneracy of the levels of the exact spectrum is an indication that the mean-field 
state spontaneously breaks the reflection symmetry of the trap: for N — + oo there are two 
degenerate condensed states in which a vortex sits either on one side or on the other side of 
the trap (in the reflection symmetry related position). The lifting of the degeneracy of these 
two states, at finite N, has been explained in terms of the rate of tunnelling of the vortex 
between these two locations [75] . 

Restoration of Rotational Symmetry 

Noticing that the GP theory breaks rotational symmetry, an improved mean-field theory 
for rotating Bose gases has been proposed and analysed in Ref.[80]. In this theory the 
GP wavefunction is projected onto an eigenstate of angular momentum. This restores the 
property that the exact groundstate for a finite N must have definite angular momentum. 
The resulting "rotating vortex cluster" state has a lower energy than the GP state, which 
may be important in situations when N is sufficiently small. 

2.3.3 Signatures of Strong Correlations 

From the first studies of rapidly rotating atomic Bose gases [21] it has been recognized that 
the groundstates include strongly correlated phases. In particular, the exact groundstate 
of N bosons in the 2D LLL at high angular momentum, L = N(N — 1)[2T], is the bosonic 
Laughlin state [81]. This state ([79]) and its properties are described in detail in §2.41 For 
now we simply note that this is a strongly correlated liquid phase of the bosons. It has a 
vanishing condensate fraction, and does not have any long-range crystalline order (in cither 
the particle or vortex density). 

Evidence for other strongly correlated states was found in exact diagonalization studies 
of small numbers of particles [82]. As the rotation frequency is swept from 57 = to the 
centrifugal limit f2 = uj±, the groundstate undergoes a series of transitions, between states 
with different values of angular momentum, L. The angular momentum L does not increase 
uniformly but undergoes steps between certain "magic angular momenta" in a sequence that 
depends on the total number of particles N. The states at these magic angular momenta 
were shown to be strongly correlated, in the sense of small condensate fraction. Furthermore, 
these states were shown to include the Moore- Read, or Pfaffian, state at L — N(N — 2)/2 (see 
§2.4|) . and a set of other uncondensed states could be accounted for in terms of condensates 
of "composite bosons" rather than of the underlying bosons [82], 

An alternative description of the rotating groundstates was introduced in Ref. 83 . This 
work proposed a generalization of the "composite fermion" construction used in the FQHE[84l 
[24] [85] to the case of rotating bosons. See §2.4l for a definition and discussion of these states. 
In this theory, one parameterizes the A-boson wavefunction in terms of a Slater determinant 
of A single-particle wavef unctions, which represents a set of A non-interacting composite 
fcrmions (|81[) . For a composite fermion state with total angular momentum Lqf the re- 
sulting boson wavefunction has angular momentum|83j L = Lcf + A (A — l)/2. Treating 
the composite fermions simply as non-interacting particles, with a Landau level spectrum, 
leads to the prediction of a set of special values of the angular momentum L* CF at which 
one expects the non-interacting composite fermions to have stable groundstates. These so- 
called "compact states" [5SJ [53] identify a set of special values of the angular momentum for 
the bosons L* = L* CF + A(A — l)/2 and provide a parameter-free trial many-boson wave- 
function at each angular momentum. The theory is very successful at describing the exact 
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Figure 8: Groundstate angular momentum as a function of rotation rate f2 for N = 9 
bosons in the 2D LLL with repulsive contact interactions. The solid line shows the results of 
exact diagonalization studies, with plateaus at the "magic angular momenta". The dashed 
line is the parameter- free prediction from the CF theory [83]. This is obtained by finding 
the groundstate angular momentum using just the set of "compact states" of the CFs. 
These compact states account for all of the magic angular momenta that appear in the 
exact spectrum; discrepancies arise only in the positions and widths of the plateaus, which 
depend very sensitively on small energy differences. 



groundstates of rotating bosons for the systems sizes that were studied (up to N — 10): the 
values L* deduced from the CF theory account for almost all of the magic angular found 
in the exact studies; the overlaps of the trial CF wavefunctions with the exact groundstates 
are very large[53]. An example of the accuracy of the CF theory is shown in Fig. [5) 

Ref. [57] has investigated the pair correlations of the exact groundstates for small numbers 
of bosons. The results show evidence for local crystalline order. These crystalline correlations 
appear even for the Laughlin state for the system sizes studied (N < 8). Since the Laughlin 
state is the exact groundstate for all N, and it is not crystalline in the thermodynamic limit, 
these crystalline correlations are interpreted as a feature of the states of small numbers of 
particles [57], 

The exact groundstates of small numbers of particles have been related to "rotating 
boson molecule" states [88 . in studies of up to N = 11 particles [89], This description is 
based on the idea of an underlying crystalline order, combined with the projection of these 
cluster states onto states of definite angular momentum, leading to angular averaging of the 
particle density. This interpretation provides an alternative viewpoint on the magic angular 
momenta. It can account for the values of the magic angular momenta in terms of the 
shell structure of ordered crystalline arrangements of the bosons, most clearly for N = 6 
and N = 9. It would be of interest to see the results of calculations of the overlaps of the 
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variational rotating boson molecule wavefunctions with the exact groundstates. 

For small numbers of particles it is perhaps moot to discuss which of a set of approximate 
theories should be used to describe the system^ Descriptions that are apparently very 
different (liquids of composite fermions, rotating vortex clusters or rotating boson molecules) 
could, in principle, all capture aspects of the groundstate and each could give useful physical 
insight. The relative merits can only be judged in their quantitative successes. However, 
qualitative distinctions between different descriptions can arise in the thermodynamic limit 
N -> oo. 



2.4 Strongly Correlated Phases 
2.4.1 The Filling Factor 

In order to make sharp statements about phases of matter, one must consider the thermo- 
dynamic limit, N — > oo. There are different ways in which to take this limit for a rotating 
gas of bosons, depending on how one chooses to scale the angular momentum L with the 
number of particles N. 

The simplest situation to consider is N — > oo keeping the number of vortices fixed. Since 
N v ~ L/N this requires L — > oo such that L/N — > const. The particle density then grows 
as n oc N, so for the mean-field interaction gn to remain finite (for example, to stay in 
the 2D LLL regime), one should take gN — > const. With this set of scalings, the Gross- 
Pitaevskii theory has been shown to be exact [73] : the groundstate is a pure condensate in 
this thermodynamic limit. We have shown explicit examples of how this limit emerges from 
the exact spectra at finite N in §2.3.21 

Although this limit is appropriate for the description of atomic BECs in many experi- 
mental situations, from the perspective of many-body physics it is not the most interesting 
limit, leading as it does to states that are fully condensed. A much more interesting limit 
can be found by exploiting analogies with the FQHE. This regime may be of importance in 
future experiments at high vortex density. 

The FQHE refers to a set of strongly correlated quantum phases of electrons in the 2D 
LLL states [531 [23]. These states are characterized by the electron filling factor 

where n D is the 2D number density of electrons. From the above mapping of the rotating 
atomic gas fjlSJt. the analogous quantity is [22] 

" n2d ^F = n2d 2MU (65) 

where the 2D density ri2d is related to the 3D density by the integral along the rotation axis, 
ri2d = J n(r)dz. From (jlOp . the filling factor can be written in terms of the vortex density 
as 

v=m. (66) 

In an inhomogcncous trap, Eqn. (|66[ ) should be interpreted as a local relation, defined on 
lengthscales larger than the mean vortex spacing. Considering the uniform limit, in which 
the particles are uniformly distributed (on average) over an area containing N v vortices, the 
relation takes the simple form 

<67 » 



5 When the exact groundstate is known (as for the Laughlin state at L = N(N — 1)), discussions reduce 
to considerations of the properties of this state. 
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Note that, for system with a large number of vortices N v , which shall be the focus of the 
studies described below, the rotation frequency is Q ^ u>± (|93|) , and so n v ~ l/(ira,j_) in 
(ED. 

It is interesting to study the nature of the groundstate of a rotating atomic gas for 
TV — » oo, -ZVv — > oo such that v — > const. A sharp question can be formulated in this 
thermodynamic limit: what is the phase diagram of a rapidly rotating atomic g as as a 
function of the filling factor vl This is the question that was first raised in Ref. [22] . and 
that we shall address in the remainder of the section. 

Note that for bosons in the 2D LLL, the filling factor can take any value, including values 
larger than one. (This contrasts with the case of fermions where the Pauli exclusion principle 
limits the occupation.) The only limitation to the theory as presented is that, for the single 
particle states to remain in the 2D LLL we require that interactions (and hence chemical 
potential) remain sufficiently small (|21[) . Since p ~ v>Vq , if Vq is small compared to the trap 
level spacings (hu)±, hu)\\) the filling factor can be very large and the system still be in the 
2D LLL. This is the situation in the experiments of Ref. [90] . where the LLL condition is 
achieved p < 2fru>± with a large filling factor v ~ 500. 

It has been argued in Refs.|221 H9, 50 that for filling factor above a critical value v c the 
groundstate is a triangular vortex lattice, so GP mean-field theory is at least qualitatively 
accurate. Indeed, as described above, the GP theory is exact in the limit v — > oo (this follows 
from the thermodynamic limit in which N — > oo with N v finite). However, for small values 
of v, corrections to mean-field theory are large. These corrections cause mean-field theory 
to fail qualitatively for v < v c , and the groundstate is replaced by strongly correlated phases 
very different from the vortex lattice [2"2"]. A very instructive way in which to understand the 
importance of the filling factor on the mean-field groundstate is to evaluate the quantum 
fluctuations of the vortices. 

Quantum Fluctuations of Vortices 

Consider the dynamics of a single vortex line in a 2D fluid (i.e. a straight vortex line). 
We shall describe the quantum dynamics of the vortex line using a canonical quantization 
approach, but the results agree with those found in other wavs |911 192) . The classical dy- 
namics of a 2D vortex, at a position X and Y in an external potential V(X, Y), follows from 
the standard Magnus force dynamics of a vortex line in a classical fluid 

-p s nY + F x xt = (68) 
+ Ps nX + F^ t = (69) 

where p s is the mass density (per unit area) of the fluid and k the circulation of the vor- 
tex. The only amendment for a quantized vortex in a superfluid is that the circulation is 
quantized, k = h/M, so we may write 

p s K = (n 2 dM)^ = n 2d h . (70) 

A Lagrangian that reproduces this classical dynamics ifH 

L = n 2 dh XY - V{X, Y) [F cxt = - W] . (71) 

Constructing the momentum conjugate to the particle co-ordinate X and applying canonical 
quantization, leads to 

QT 

H x = =n 2d hY (72) 

oX 

[X,fl x ]=ih => [X,Y] = —^—. (73) 

27rn2d 



3 This is easily checked by constructing the Euler-Lagrange equations for J7TJ. 
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The X and Y co-ordinates are conjugate^ so obey the generalized uncertainty relation 



AX AY > (74) 

- 47m 2d 

which implies 

AX 2 + AY 2 > — !— . (75) 
27rn 2 d 



The result ([75]) makes physical sense: one cannot locate the vortex line to a distance 
less than the mean 2D separation between the particles. It is interesting to note that this 
result has an entirely classically interpretation, but emerges from a quantum calculation 
due to the cancellation of Planck's constant in the circulation with Planck's constant in the 
commutator. Furthermore, the result (|75[) is consistent with the calculation of fluctuations of 
the vortex (f60|) for the exact one- vortex wavefunction (f59|) . Eqn. ([60]) implies AX 2 + AY 2 ~ 
£ 2 /N ~ l/w 2 d, where we take the typical particle density n 2 d ~ N/£ 2 noting that the N 
particles are within an area of ~ I 2 in this inhomogeneous state (|59p. 

The importance of the filling factor for the properties of a large vortex lattice becomes 
clear if one applies a form of Lindemann criterion and asserts that the vortex lattice will 
become unstable if the (rms) quantum fluctuations in vortex position are larger than some 
multiple qj, of the vortex spacing 



^AX 2 + AY 2 = * > a L xa v = a L i -2— (76) 
V27rn 2 d y V3n v 

_ n 2 d , \/3 , , 

v = > v c = T . (77) 

Putting in a typical value for the Lindemann parameter, a 2 L ~ 0.02 93 , one finds v c ~ 7. 

The theory described above for a single vortex neglects the kinetic energy of the vortex. 
It amounts to the evaluation of the quantum fluctuations of only the guiding centre co- 
ordinate. Including the inertia of the vortex leads, in addition, to a zero-point cyclotron 
motion of the vortex position. This zero point cyclotron motion gives an equal contribution 
to the quantum fluctuations of position as that of the guiding centre. Combining these 
fluctuations, one finds 

AX 2 + AY 2 > — !— (78) 
7rn-2d 

that is, twice the value ([75]) from the guiding centre fluctuations alone. Using this expression 
([78]) within the Lindemann analysis leads to a critical filling factor that is twice as large, 
v c = V3/ (2na 2 L ), which evaluates to v c ~ 14 using the same value for the Lindemann 
parameter. This is the estimate that was given in Refs. [9"31 12"2"] . 

A more controlled calculation of the quantum fluctuations of a vortex lattice involves 
the consideration of the collective modes of the entire array, rather than just the motion 
of a single vortex. This calculation has been performed for a uniform vortex lattice ar- 
ray, by determining the low frequency "Tkachenko" modes [94] [95] [96] and evaluating their 
contributions to the quantum fluctuations of an individual vortex [IS]. Applying the same 
Lindemann criterion as above, leads to a critical filling factor of v c ~ 8. A more complete 
calculation that includes, in addition to the Tkachenko modes, the gapped "inertial modes" 
of the vortex lattice[50j leads to a critical filling factor of v c ~ 17. The difference between 
these two results, of roughly a factor of two, is again due to the neglect or inclusion of the 
inertial contributions to the fluctuations of the vortices [50]. That the results of calculations 
based on the collective modes are so close to those found from single-vortex calculations is 
an indication that the fluctuations are dominated by the short-wavelength modes. 



7 These are the guiding-centre co-ordinates of a particle in a single Landau level. 
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Figure 9: Numerical evidence for the development of short-range crystalline order at filling 
factor v = N/N v > 6 [52], for bosons in the 2D LLL on the torus geometry with aspect 
ratio chosen to be commensurate with the triangular vortex lattice. The graph shows the 
excitation energies at a set of different points (K Xl K y ) in the Brillouin zone as a function 
filling factor v. For v > 6 the excitations at reciprocal lattice vectors of the triangular vortex 
lattice (filled symbols) become very small. This is a signal, in the finite-size system, of a 
tendency to translational symmetry breaking. (For v < 6 there is a low energy excitation 
which is located at a position in the Brillouin zone that is consistent both with the Read- 
Rezayi states and with the "smectic" phase.) [Reprinted figure from: N.R. Cooper, N.K. 
Wilkin, and J.M.F. Gunn, Phys. Rev. Lett 87, 120405 (2001). Copyright (2001) by the 
American Physical Society. ] 

2.4.2 Numerical Evidence for Crystalline Phases 

While instructive, considerations based on the Lindemann criterion are hardly predictive, 
depending very sensitively on which, itself, is estimated [93j from the thermal melting of 
3D crystals! For this reason it is useful to have a direct determination of the transition. 

The transition to the vortex lattice phase at large v has been studied in large scale 
exact diagonalization studies 22, 97J. The strategy is to work on a geometry with periodic 
boundary conditions (the torus geometry, §B.3|) . which is consistent with the formation of 
a vortex lattice. The signal of crystallization is the collapse to very low energies (above the 
groundstate energy) of a set of excitations at momenta that are reciprocal lattice vectors of 
the vortex lattice. (This is directly analogous to the signal of rotational symmetry breaking 
in a trap, leading to quasi-degeneracies of energy levels at angular momenta spaced by 
AL, which was discussed in ^2.3.21 ) By looking for the emergence of broken translational 
symmetry, it was estimated that the transition to the triangular vortex lattice occurs at 
v c ~ 6 [22 . Numerical evidence for the appearance of a triangular vortex lattice is shown 
in Fig. At v = 6 the largest system that could be studied by exact diagonalization was 
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Figure 10: Schematic diagrams of the "smectic" phase that appears in exact diagonalization 
studies on the torus (left), and its possible effect in disordering the vortex lattice in a trapped 
gas (right). The vortex lattice is quantum-disordered along lines of the vortices, giving a 
phase with broken translational symmetry in just one direction [97). 



TV = 48, N v = 8. Since the lattice contains only 8 vortices, there are surely significant finite 
size effects in the estimate of the location of the phase transition. 

More recent work [97] focusing on v = 2 has shed more light on this issue. In the 
regime v — 2 — > 6 the numerical results were shown to be consistent with the existence 
of a competing crystalline phase with "smectic" order - that is, a stripe state with broken 
translational order in only one direction. This can be viewed as a phase in which the vortex 
lattice is quantum melted due to fluctuations of vortices along lines, leading to averaged 
density patterns of the form shown in Fig. [TUJ While this state appears to describe the 
groundstate in the system sizes that can be studied numerically (N v < 12), it is possible 
that finite-size effects may favour this over competing phases. Estimates suggest that the 
finite-size effects in experiments on cold gases will mean that the smectic phases are likely to 
be important in practice[97J. The first indications of quantum disordering of vortex lattices 
as filling factor is reduced is likely to be the appearance of (local) stripe ordering, perhaps 
as in Fig. CEl 

In the range of filling factor 2 < v < 6 close to the appearance of the short-range 
vortex-lattice ordering at v ;S 6, the usefulness of exact diagonalization studies becomes 
limited with respect to determining the groundstates in the thermodynamic limit, owing 
to the influence of finite-size effects. This is apparent from the strong dependence of the 
numerical spectrum on the boundary conditions (aspect ratio of the torus), which shows 
there to be competition between different phases: the vortex lattice, Read-Rezayi states 
and the smectic state. Understanding the nature of the groundstates in this regime is a 
very difficult theoretical challenge, involving a strongly interacting many-body system close 
to a quantum phase transition. This would be a particularly interesting regime to explore 
experimentally. 

2.4.3 Incompressible Liquid Phases 

For v < v c the vortex lattice phase (a triangular vortex lattice) is unstable to quantum 
fluctuations, and is replaced by a series of strongly correlated phases. These phases are best 
understood at small filling factors, far from the transition to the vortex lattice phase. As 
the filling factor approaches this transition, our understanding becomes much poorer (see 
the discussion in 8 32. 4. 2|) . 
Laughlin State 

The strongly correlated phase that is best understood is the bosonic Laughlin state [51], 
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at v = 1/2. The wavefunction describing this state is 



N 

*Laughlin({C*})«n^-0) 2 - ( ?9 ) 

i<j 

This is the exact groundstate for contact repulsion (|13[) at total angular momentum L = 
N(N— 1)121]. This may be noted from the fact that this is the unique symmetric polynomial 
function of N variables that has the property that it vanishes when any two co-ordinates 
coincide - it is therefore the unique zero energy eigenstate of the contact interaction Q13p in 
the 2D LLL. 

The state has the property that its average density is uniform [81], up to a radius ~ 
2\J~Nl beyond which the density falls to zero. Although the particle density of this state is 
uniform, there are still vortices in the system. However, unlike the vortex lattice phase ([40]) . 
these vortices are not localized in space (translational symmetry is not broken). Rather, 
the vortices can be viewed as being bound to the particles[98 , 99 . This can be seen by 
noting that the wavefunction (|79|) changes phase by 2 x 2ir each time the position Q of any 
particle i encircles the position of any other particle Q=^i ■ Each particle therefore experiences 
two vortices bound to the position of every other particle, so the total number of vortices 
experienced by any one particle is N v = 2(N — 1). Thus, (|T9"|) describes a phase that has 
filling factor v — N/N v — 1/2 in the thermodynamic limitd 

The Laughlin state for bosons shows all of the characteristics familiar from the conven- 
tional FQHE[23l [24]: it is an incompressible fluid, with gapped collective excitations in the 
bulk, and gapless edge modes; the particle- like excitations have fractional particle number 
and fractional statistics (they are "abelian anyons" ) [23] . 

Incompressibility of the bosonic Laughlin state for rotating bosons is evident from exact 
diagonalization studies of the collective excitation spectrum. The spectrum shows the char- 
acteristics of an incompressible liquid (see Fig. II 1[) in which the excitations are gapped at all 
momenta rather than having the linear sound mode of a compressible liquid. Calculations 
of the energy gap on finite size systems, show convergence to the thermodynamic limit, with 
a collective mode gap of 0.095(5) x (27rV ) [1001 DUE] E 

For L > N(N — 1), the groundstate of the contact interaction becomes degenerate, with 
a subspace spanned by quasi- hole excitations of ([79]) [102] . These states may be written 

N Wqh 

Kl) OC ]J II (6 - W «) X ^Laughlin , (80) 
i— 1 a — 1 

where w a are the co-ordinates of _/V q h quasiholes0 From general properties of the Laughlin 
states [23], these quasiholes have a fractional particle number, giving rise to a density deple- 
tion the integral of which is exactly 1/2 of an atom; they behave under adiabatic exchange as 
particles with fractional exchange statistics |103j ("anyons") with statistics parameter 1/2, 
corresponding to "semion" statistics. 

Incompressibility implies that the chemical potential has a discontinuity: the chemical 
potential for adding particles to a bulk system, is not the same as the chemical potential 
for removing particles, Since the chemical potential for a bulk system is the derivative of 
the energy density with respect to the number density, this implies that the energy density 



8 The exact relation N v = 2N — 2 should be interpreted as N v = N/v — S where S is the "shift" of the 
state on the sphere, see i]B.2l 

9 The numerical results in Refs. [T00l I101| are stated to be in units of \fly2Trhu) ±a s / 'ay = 47rVb- For 
consistency with the numerical calculations of Ref. 22 i and here, I assume that the results of Refs. [T00lll01| 
are, in fact, in units of 2ttVq. 

10 As written the states H80H are not eigenstates of total angular momentum. They form an over-complete 
basis for the zero-energy states with L > N(N — 1). 
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Figure 11: Excitation spectrum of the Laughlin state at v = 1/2, for bosons interacting 
with contact interactions, obtained by exact diagonalization studies for TV = 8 and N v = 16 
on a torus with square geometry, a/b — 1. The states are labelled by the magnitude of 
the conserved wavevector in units of 2ir/a (see TO.3|) . The zero energy state at \K\ = is 
the Laughlin groundstate; the excitations are gapped at all wavevectors, characteristic of an 
incompressible liquid. 



has a cusp for an incompressible state. One can relate these thermodynamic quantities 
to the creation energies of quasiparticles, A qp , and quasiholes, A q h, in the incompressible 
state. For the v = 1/2 Laughlin state, quasiparticle and quasihole excitations have fractional 
particle number of 1/2. Thus, adding a single particle to the bulk of the system involves the 
creation of two quasiparticles, and an energy increase of /i + = 2A qp ; removing a particle 
involves the creation of two quasiholes, and an energy increase of — /x_ = 2A q h- Hence, the 
discontinuity in chemical potential is /i + — /i_ = 2(A qp + A q h). For the Laughlin state with 
contact interactions, the energy to create a quasihole vanishes, A q h = 0. The quasiparticle 
creation energy is non-zero, and is available only from numerical studies 22, 100, 101 . The 
numerical results are consistent with the expectation that the quasiparticle energy is equal 
to the collective mode gap jl00[ 1 101] . [The large wavevector collective mode can be viewed 
as a widely separated quasiparticle quasihole pair, so is a measure of (A qp + A q h). From 

FignUthiS iS A qp + A qh = A qp ~ 0.1(271%)-] 

Composite Fermion States 

At certain higher filling factors, v > 1/2, there appear strongly correlated phases that are 
accurately described in terms of non-interacting "composite fermions" [83j I100[ 11011 1104] . As 
described in ij2.3.31 the states formed from non-interacting composite fermions provide an 
accurate description of the groundstates of small numbers of rotating bosons in a trap. These 
are the small-system signatures of bulk incompressible liquid phases (in the thermodynamic 
limit). 

Composite fermions are formed from binding vortices to particles. Since the underlying 
particles are bosons, one can form a composite particle which has Fermi statistics by binding 
a single vortex to the location of each particle. This may be made explicit within the Jain 
construction [84], by writing the many-particle wavefunction for the bosons as|83j 

*({r fe }) = PlllII^ " 0) ^CF({r fc }) • (81) 

i<j 

The Jastrow factor rXi<j(Ci — Cj) causes any particle i to experience a single vortex at 
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the location of any other particle j 7^ i. Since this factor is completely antisymmetric 
under particle exchange, to obtain a bosonic wavefunction the function ^cf must also be 
antisymmetric: this is the wavefunction of the composite fermions. In general terms, one can 
appreciate why this is a useful variational state: the Jastrow factor suppresses the amplitude 
for two particles to approach each other, as does the composite fermion wavefunction. Taken 
alone, these two factors would give a bosonic wavefunction that is a zero energy eigenstatc 
of the contact interactions. However, for L < N(N — 1), this construction requires ipcF 
to include basis states that are not in the lowest Landau level. To recover a bosonic state 
within the LLL, one must project the single particle states into these states, as represented 
by the operator Plll[H]- For L < N(N — 1) the projected wavefunction does not vanish as 
two particles approach each other, so it is a state with a non-zero contact interaction energy. 

Since one vortex is bound to each particle, in a large (uniform) system the composite 
fermions experience 

n^ F = n v - n 2d (82) 
vortices per unit area, so have filling factor (f65|) 

CF _ n 2d 



v 



(83) 



Treating the composite fermions as non-interacting particles which completely fill p Landau 
levels, i.e u CF = ±p, one is led to the bosonic version of the "Jain" sequence 

■ " (84) 



p± 1 



At these filling factors, the composite fermion theory predicts the appearance of an incom- 
pressible liquid state of the bosons, and provides a trial wavefunction (|8ip . 

The states constructed in this way have been shown to successfully describe the results 
of exact diagonalization studies in a variety of ways. 

For trapped systems (in the disk geometry): (i) The prediction of the "magic angular 
momentum" for states of small numbers of rotating clusters [53], (ii) Large overlaps of the 
groundstate wavefunctions with the composite fermion wavcfunctions [55] . (iii) The form of 
the edge excitation spectrum for a bulk incompressible liquid of this kind is known 1051 1106] . 
The edge modes found in exact diagonalization studies are consistent with these expectations 
for bulk regions at v = 1/2. 2/3. 3/4(107]. 

For uniform systems, the evidence of composite fermion states at v = 1/2,2/3 and 
3/4 includes: (i) The existence of a sequence of states for different particle number N 
at the expected "shift" on the sphere [see Eqn. (|132|) ]. which have a uniform groundstate 
charge density and an energy gap (i.e. consistent with incompressible liquids) [100, 101 . 
(ii) The prediction of the quantum numbers of the low- lying excited states [101] ■ (iii) Large 
overlaps of the exact groundstate and of the low-lying excited states with the parameter-free 
wavefunctions formed from the above composite fermion wavefunctions [T04] . 

Extrapolation of the energy gap at v — 2/3 shows a value of about 0.05 x (2-7rVo) in the 
thermodynamic limit [101] F*l a similar extrapolation for v = 3/4 is unreliable owing to the 
small number of system sizes available. 

Moore-Read and Read-Rezayi States 

One of the most interesting aspects of the physics of rapidly rotating Bose gases is the 
prediction [51] [551 [108, 109, 97 of the appearance of non-abelian phases: incompressible 
liquid phases whose quasiparticle excitations obey "non-abelian exchange statistics" . 

Non-abelian phases of matter are currently attracting intense theoretical and experi- 
mental interest in the condensed matter physics community [110] . In part, this is due to 



'See footnote 9 on page I29l regarding the units in Ref. 1001 1101] . 
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the possibility to explore experimentally this most exotic consequence of quantum many- 
body theory. In part, these studies are motivated by the possibility ultimately to build a 
"topological quantum computer" , which makes use of the properties of a non-abelian phase 
to form a quantum register with topological protection against decoherence. Experimental 
work is largely centered on studies of the v = 5/2 FQH state of electrons in semiconductors. 
This state is thought likely to be a non-abelian phase of matter described by the Moore- 
Read, or Pfafhan, state [111] or a close relative thereof. Experiments are working towards 
the detection of the expected unconventional exchange statistics of the quasiparticle excita- 
tions of this state. Unfortunately, the creation of a universal topological quantum computer 
would not be possible with the Moore-Read state, as its properties do now allow arbitrary 
unitary transformations of the quantum register. Other non-abelian phases, known as the 
Read-Rezayi states [TT2L would allow universal topological quantum computation. There is 
some theoretical evidence that these states might appear in semiconductor systems, but this 
is still very preliminary [113]. 

With this background, it is striking that theoretical studies have shown very convincing 
evidence for the appearance of the Moore-Read state [111] and the Read-Rezayi states [112] 
for rotating Bose gases with realistic two-body interactions. One of the major goals of the 
creation of rapidly rotating atomic gases would be to allow (possibly the first) measurements 
of non-abelian phases of matter. 

The construction of Moore-Read and Read-Rezayi states for bosons may be viewed as a 
generalization of the bosonic Laughlin state. The Laughlin state for bosons is the densest 
exact zero energy eigenstate of the 2-body contact interaction (|13p within the 2D LLL. Sim- 
ilarly, the Moore- Read and Read-Rezayi states are the densest exact zero energy eigenstates 
of a fc + 1-body contact interaction 

N 

S(r il -r i2 )S(r i2 -r i3 )...S(r ik -r ih+1 ). (85) 

»1<»2<".«* + 1 = 1 

For N divisible by fc, the groundstate wavefunctions of (|85p may be written in a simple way 



N/k N/k 

n cc< - o) 2 n - c 

i<jeA KmeB 



(86) 



where S symmetrizes over all possible ways of dividing the N particles into the fc groups 
(A,B...) of N/k particles each|114. It is straightforward to convince oneself that ([56"]) 
vanishes when the positions of fc + 1 particles coincide, as required. By counting the degree 
of the polynomial L — N(N/k — 1), or the number of vortices N v — 2(N/k — 1) experienced 
by each particle, one sees that these states describe bosons at filling factor 

^ = ~ ■ (87) 

fc = 1 is the Laughlin state; fc = 2 is the Moore-Read state; fc > 3 are the Read-Rezayi 
states. Although these states have a similar structure, it should be remembered that the 
physics of these phases is very different - with quasiparticle excitations that obey abelian 
(fc = 1) or non-abelian (fc > 2) exchange statistics. 

Numerical evidence for the Moore-Read state (fc = 2) for bosons interacting with contact 
interactions has been reported on the disc |115j . torus [2"2] and spherical [1 001 1 1 04j geometries. 
Large overlaps of the exact wavefunctions with the model states ([86]) are found. While a 
large wavefunction overlap is encouraging, this is not necessarily the best way to characterize 
a phase of matter. In the thermodynamic limit, owing to the exponential increase in the size 
of the Hilbert space, the wavefunction overlap with any trial state will surely vanish, even 



32 




Figure 12: Low energy spectrum for the v = 1 state with contact interactions, obtained by 
exact diagonalization studies on a torus with square geometry for N = 16 and iV v = 16. The 
states are labelled by the magnitude of the conserved wavevector on this periodic geometry. 
The appearance of low-energy states (shown as filled circles) at K = (8,8), K = (8, 0) and 
K = (0, 8) (the last two are indistinguishable in this plot by symmetry) is consistent with 
the expected three-fold degeneracy of the Moore-Read state. [K is measured in units of 
2-k/o, where a is the side of the (square) torus.] 



if the two wavefunctions describe the same topological phase. A robust characterization 
of the topological phase is provided by the groundstate degeneracy on the torus [TIB"] ; this 
degeneracy is expected to survive (and even improve) in the thermodynamic limit, provided 
the wavefunction is in the same topological phase as the trial state (j56"|). For the Moore- 
Read state, this degeneracy appears clearly in the spectrum for the system sizes that can 
be studied numerically [2"2l [TT7] . see Fig. [121 

Evidence for the appearance of phases of matter that are described by the Read-Rezayi 
states (k > 3) at v = k/2 for contact repulsion (| 13[) has been found on the torus geometry [22], 
While wavefunction overlaps are large, at least for k that is not too large, the expected 
groundstate degeneracy is less clearly resolved than for the Moore-Read state (k = 2) 108 . 
Studies of the groundstates in the spherical geometry at v = 3/2,2, at the shifts [see 
Eqn. (|132p ] appropriate for the k — 3 and k — 4 Read-Rezayi states, were determined to 
be inconclusive [ 1 OOl 1118] with no clear sign that the thermodynamic limit had been reached. 
It thus appears that, for contact interactions, the correlation lengths of the groundstates 
at v = 3/2 and v = 2 are somewhat larger than the available system sizes. (The ori- 
gin of this lack of convergence is related to the competing "smectic" state at intermediate 
filling factors[H3) N2.4.21 ) However, it has been found that by changing the inter-particle 
interactions to introduce an additional small non-local repulsion, the spectra at v = 3/2 
and 2 smoothly evolve into spectra with very clear groundstate degeneracies which show 
excellent evidence of being in phases described by the k = 3 and k = 4 Read-Rezayi 
states [1Q81 11091 197] . It appears that the small non-local interaction makes the correlation 
lengths smaller than the system sizes available. These results are described further in ^2.4.51 

2.4.4 Feshbach Resonance 

For a model of rotating bosons interacting by a Feshbach resonance, it can be shown that (for 
certain parameters) the exact groundstates are the Moore- Read and Read-Rezayi states [1191 
[8] . The hybridization of atoms into molecules leads to a suppression of the repulsive inter- 
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atomic interactions. Tuning the resonance to the point where the net two-body repulsion 
vanishes leads to an effective theory for dressed atoms (which are resonating into molecules) 
in the 2D LLL which involves effective three- and four-body contact interactions. These 
arise from processes in which two atoms hybridize into a molecule which then collides with 
an atom, or another molecule. 

2.4.5 Dipolar Interactions 

The effects of dipolar interactions on the strongly correlated phases of rotating bosons have 
been investigated in Refs.[62l H08|. l97l I120|, I117j . Qualitatively similar effects are found 
from studies of a model interaction in which a non-zero V2 Haldane pseudo-potential is 
introduced [1211 [1091 H20| in addition to the contact Vq. In this section we will focus on the 
case of dipolar interactions, with pseudo-potentials (|4*l)j) . 

The strength of the dipolar interaction relative to the contact interaction is characterized 
by the ratio a = V2/V0 (|49j) . One expects the nature of the groundstates to vary with a. 

For small positive values of a the nature of the groundstate does not change significantly. 
However, it has been shown that this can lead to improvements in the overlaps of the 
groundstates (as computed in exact diagonalization studies) with the Moore-Read state at 
v = 1 and the Read-Rezayi states at v = 3/2 and 2[TUH1 Ml HH EH]- lt appears that 
the effect of a small amount of additional long-range repulsion is to reduce the correlation 
length of the groundstate sufficiently that it is less than the available system sizes. Numerical 
results then show convergence to the thermodynamic limit. Evidence for the k = 3 Read- 
Rezayi state at v = 3/2 is shown in, Fig. [13l Currently, these results provide by far the 
most convincing theoretical evidence for the existence of a Read-Rezayi topological phase 
in a system with realistic two-body interactions. 

For larger values of the non-local interaction, the nature of the groundstate changes 
character. 

The change in the nature of the groundstate at v = 1/2 was described in Ref. [62] . 
based on the results of exact diagonalization studies on the torus. For a = the exact 
groundstate is the Laughlin state. For a Z 0.5 the groundstate changes character abruptly 
(the overlap with the Laughlin state becomes very small). The states for a <; 0.5 were found 
to be compressible crystalline phases. Numerical evidence was given for a stripe phase at 
a = 0.528, and for the q = 4 bubble crystal phase at a = 0.758. These results are consistent 
with expectations of mean-field theory, Fig. [H which indicates stripe and bubble crystal 
phases at these values of a. Mean-field theory is expected to be reliable only for sufficiently 
large v. Why does it work so well even for v = 1/2 in this regime of a > 0.5? The answer 
can be found by noting that the mean-field phases at these values of a involve clustering of 
particles, so the quantum fluctuations are reduced. For example, the q = 4 bubble crystal 
phase at a = 0.758 has vq = 2 particles per bubble, so can be viewed as a crystal of pairs of 
particles. A calculation of the Lindemann criterion for quantum melting of the q — 4 bubble 
crystal phase leads to a critical filling factor of v^ =A ~ 0.4, using a Lindemann criterion that 
gives v c ~ 6 for the critical filling factor at which the triangular vortex lattice is unstable 
to quantum fluctuations [62 . This estimate is consistent with the observation of crystalline 
order even down to v = 1/2. 

The collapse of the incompressible states for large values of a has been reported at filling 
factors v = 1 [120[ 1117] . v = 3/2 |108j and v = 2[97]. In all cases the transition is into 
compressible states. At v = 1 the Moore-Read state is destroyed for a £ 0.4, and evidence 
has been given for the appearance first of a stripe phase [12"0|,lll7] and then of bubble crystal 
phases [120]. These results are consistent with the expected increased stability of the stripe 
and bubble crystal phases at v = 1/2 as the filling factor is increased 108J. (The q = 4 
bubble crystal phase is a crystal of clusters of vq = 4, 6 and 8 particles for v = 1,3/2 and 2 
respectively. ) 
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Figure 13: Numerical evidence for the appearance of the k = 3 Read-Rezayi state for 
bosons at v = 3/2. (N = 18, 7V V = 12 on a torus with two different units cells, for contact 
interactions with a small amount of dipolar interaction, a — V2/V0 — 0.380.) The spectrum 
depends weakly on the unit cell, so the results are representative of the thermodynamic limit. 
There is a very clear 2-fold degeneracy of the groundstate, exactly as expected for the k = 3 
Read-Rezayi state. (The overall 4-fold degeneracy is recovered by the additional two-fold 
centre-of-mass degeneracy due to the half-integer filling factor.) From Ref. |108j . [Reprinted 
figure from: E.H. Rezayi, N. Read, and N.R. Cooper, Phys. Rev. Lett 95, 160404 (2005). 
Copyright (2005) by the American Physical Society. ] 
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As compared to the case of contact interactions, one qualitatively new effect of small 
non-zero value of a is that it allows the emergence of other strongly correlated phases 
for v < 1/2. At v = 1/3 numerical evidence has been given [TIT] for the appearance of 
the composite fermion liquid [122] state. One further anticipates the appearance of other 
incompressible states (e.g. the Laughlin states at v = 1/m with m > 4) and the transition, 
at sufficiently small v, to a Wigner crystal phase of the individual atoms, analogous to 
that for electrons in FQHE systems or dipolar interacting fermions[123j. The properties of 
rotating dipolar bosons in this regime remain to be studied in detail. 



2.5 Density Distribution in a Trap 

In much of the above, we have treated the rapidly rotating gas in the uniform limit. This 
applies when £1 — uj± and the residual transverse trapping potential in the rotating frame 

V i i(r ± ) = hf(u,l-n 2 )rl (88) 



vanishes. In this formula r± = \J x 2 + y 2 is the radial distance from the rotation axis. 
Clearly £1 = oj± requires the numbers of vortices N v and of particles N to be infinite. In 
any finite system one has f2 < lu± so there will remain a weak harmonic trapping potential. 
This trapping leads to inhomogeneity of the density, which we now discuss. 

The first corrections to the particle density at large but finite N v can be found by using 
the local density approximation (LDA) . This amounts to the assumption that at each point 
in space there is a locally defined chemical potential 

IM(r±)=fi-V£(r ± ). (89) 

Then, the 2D particle density is found from the local chemical potential, as 

n 2 d(rx) = n 2 d[fJ.(r±)} , (90) 

where ri2d[/-t] is the 2D density of a uniform system at chemical potential fi. Within a 
reasonable assumption of locality in the physics determining the equation of state, one 
expects the LDA to be accurate provided the particle density varies slowly on a lengthscale 
of the correlation length of the relevant phases; for the situations we have been discussing 
this correlation length is of order the mean vortex spacing a v ~ a±. 

Through (|90|) . the spatial density profile n2d(f±) depends on the equation of state, 7i2d[/-t]- 
The profile therefore takes different forms in the vortex-lattice (y > v c ) and strongly corre- 
lated (y < u c ) regimes. 

A useful general result exists. For a Bose gas with repulsive interactions, the density 
falls to zero for fx < 0. Therefore, the edge of the cloud is at a radius 



where /i(0) is the chemical potential at the centre of the cloud. Since the average vortex 
density is n v — 2MQ/h, the total number of vortices within the Thomas-Fermi radius R± 
is 

o 2irMn 2u(0) 

Inverting this relation, and taking the limit A^ v 3> 1 in which the LDA is valid, shows that 
the rotation frequency is 

SL * + 0{l/Nl). (93) 

Cl>_|_ Tlluj_ iv v 
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Thus, for 7V V sufficiently large (larger than Utu±/ fi(0)), the rotation frequency is close to the 
trap frequency u>±. 

Deep in the vortex lattice phase v ^> 1, the equation of state can be found from mean- 
field theory, Eqn. ([45]) . giving /i = 2(3avVq where v is the local filling factor ([66]) . The local 
density approximation therefore predicts a density distribution of the form [43] 



This is the inverted parabola familiar from the Thomas-Fermi approximation for a trapped 
BEC. This been shown to provide an excellent description of vortex lattice phase (v ^> 1) 
within the 2D LLL[?31 [UJ, showing that the LDA is accurate in this regime. Further 
applications of the LDA allow one to relate /z(0) to the total number of particles N, and fl 
to the angular momentum 2TJ. This leads to the result that the number of vortices 
inside the Thomas-Fermi radius scales as 



for large N v . Detailed theoretical studies of the particle density and vortex distributions of 
finite vortex lattices in the mean-field LLL regime have been reported in Refs.|43, 4lll42 l[T24] . 
Among the motivations for these considerations was the observation that if the vortices 
adopt a uniform triangular array the mean particle density follows a Gaussian profile on 
average[38], rather than the more conventional Thomas- Fermi profile. The reconciliation of 
these views led to the conclusion that the vortex lattice distorts slightly towards the edge 
of the trapj43][42][4T], in such a way that the above Thomas- Fermi profile is achieved. The 
arrangement of vortices outside the Thomas-Fermi radius has been discussed in Refill 24 . 

In the strongly correlated regime, the equation of state has a jump in the chemical 
potential at the density of each incompressible phase. As a consequence, within the LDA the 
filling factor does not fall smoothly with radius. The density distribution shows a "wedding 
cake" structure [125] characteristic of incompressible phases, with regions in which the 2D 
particle density is pinned to values v*n v where v* arc filling factors at which the groundstate 
is incompressible. This is expected to give clear experimental evidence for the appearance 
of incompressible phases in images of the density distribution following expansion|125 . 

2.6 Quantum Melting in the 3D LLL Regime 

The 3D LLL regime (|35l) can be relevant for a rapidly rotating Bose gas in an anisotropic trap 
Ci/ii -c u;_i_[38]. The regime is also of importance for the situation [126] in which the rotating 
gas is sliced into many parallel layers by a ID optical lattice with wavcvcctor parallel to the 
rotation axis, see Fig. [14] The optical lattice leads to an increased effective mass My for 
motion along the rotation axis, My > M. For a deep lattice in the tight binding limit, with 
lattice period d\\ and single particle tunneling amplitude J, the effective mass for motion 
along the rotation axis is 



which diverges as J — ► 0. Hence, the oscillator energy wii oc I/Mm is reduced, leading to an 
anisotropic confinement cJm -C luj_. 

The use of a ID optical lattice has been suggested as a means of increasing the quantum 
fluctuations of the vortices and allowing entry into the strongly correlated regime [121)1 11251 
11271 1128] . This also has the advantage of increasing the interaction energy scale of the 
strongly correlated phases (see W. 2p . 
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Figure 14: Schematic diagram illustrating the use of an optical lattice, directed along the 
rotation axis, to slice the rotating condensate into many parallel layers. Each layer has a 
reduced filling factor, z/i ayor ~ v/N^, where A^l is the number of layers. In addition, the 
tight axial confinement of each layer (ajj ayel ?S ) leads to an increased energy scale for 
interactions (f5T|) . 



The origin of the enhanced quantum fluctuations of the vortices is simple to understand 
in the limit of a very deep lattice in which the layers become decoupled. Then, the relevant 
quantity controlling the quantum fluctuations is the filling factor per layer 

N v 

^ * A^l = Al (97) 

where ~ a ||/^|| is the number of layers that the optical lattice forms, which could be 
of order Mayor ~ 50 — 100. The vortex lattice in each layer will disorder due to quantum 
fluctuations for Z4 ayer < v c (where v c ~ 6 is the 2D result [221 [49J, [50] ) , leading to quantum 
melting at a very high total filling factor v < v c Nt,. 

The quantum fluctuations of a uniform vortex lattice (N v large) in the 3D LLL regime 
of a layered geometry have been investigated theoretically in Ref. [125j . The critical filling 
factor for a uniform system was estimated as a function of the coupling between the layers 
both by applying the Lindemann criterion based on the quantum fluctuations of the vortices, 
and by evaluating the condensate depletion. The degree of interlayer coupling is controlled 
by the dimensionless ratio of the interlayer tunneling energy J and the mean interaction 
energy in a single layer, Vbfiaycr with V® ([3"Tj) evaluated for the subband thickness a[p yer < d\\ 
of each individual layer. The results show the crossover from the regime of strongly coupled 
layers J 3> Vb^layer where the 2D LLL regime is valid and quantum melting occurs for 
v ;$ u c , to the regime of decoupled layers for J <C Vb^iayor and quantum melting requires 
only v < v c N h . 

The criterion for quantum melting in the intermediate 3D regime can be understood 
in simple geometric terms. In this 3D geometry, the vortices do not simply fluctuate as 
straight lines. Quantum fluctuations involve also bending of the vortex lines - the Kelvin 
modes of the vortices [129] . The relevant Kelvin modes are those with wavelengths longer 
than a certain minimum wavelength, which one can think of as the lengthscale over which 
the vortices remain straight. This minimum wavelength is set by the healing length ([9]) 
along the rotation axis [125j 



or, if the healing length becomes less than the lattice period £ii < dii , then one should use 
di\ as the minimum wavelength. Note that, for a system in the 3D LLL ([35)) . the healing 
length satisfies £y <C aii , so the vortex fluctuations do involve flcxural motion. 

Applying a Lindemann criterion |125j to these fluctuations of the vortex lines leads to the 
conclusion that quantum melting of the 3D vortex lattice occurs when the filling factor in 
a slab of thickness £m is less than a critical value of order 1. We take this critical value to 
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be v c ~ 6 for consistency with the 2D studies [22, 49, 50 . Then the condition for quantum 
melting of the 3D vortex lattice is 



< v c , (99) 



where n is the 3D particle density. With this criterion in mind, and using (|98p , the condition 
for validity of mean-field theory in the 3D LLL regime (135)) is 



a. 



n > a s n* ~ -f . (100) 

It is only in this limit IjlOOp that one can model the 3D LLL regime by a Thomas-Fermi 
profile for the particle density (|127p along the rotation axis. 

The quantum fluctuations of a confined vortex lattice (N v finite) in a layered geometry 
have been studied in Ref. |127| I128j . This work provides a detailed study of the Tkachcnko 
modes for the finite-size vortex lattice. The quantum melting of clusters of vortices is found 
to be inhomogeneous, due to the inhomogeneous particle density. Melting can lead to a 
fluid of vortices, or to a shell structure in which shells of vortices can decouple from their 
neighbours but maintain inter-vortex correlations within a shell. A critical comparison is 
made with the expected location of the vortex-lattice/ vortex liquid boundary in the inho- 
mogeneous system as predicted by the local density approximation. For small vortex arrays 
(N v small) in the limit of weak interlayer coupling, one expects a transition to quantum Hall 
"Mott" phases in which number fluctuations between layers are suppressed |125j . 



3 Atomic Bose Gases in Rotating Optical Lattices 

The effects of a periodic pinning array on vortices in superconducting materials is a matter of 
significant experimental and theoretical interest [130[ I131L 11321 [133J. At low temperatures, 
the vortex lattice can undergo a transition from the triangular Abrikosov lattice (in the 
absence of pinning) into phases that are pinned to the lattice, with a series of different 
ordered structures. The nature of the groundstate involves the interplay between the vortex 
density n v and the density of periodic pinning sites n s it ei and can involve rich and interesting 
commensurability effects [13Tj as a function of the ratio 

f = ^-, (101) 

^sitc 

the mean number of vortices per pinning site. 

The interplay between vortex-vortex interactions and a vortex pinning potential can 
be conveniently studied in rotating atomic BECs. Regular arrays of pinning sites can be 
imposed by a rotating optical lattice potential [134j . and allow detailed investigations of 
the vortex lattice structure and the dynamics. Furthermore, the extension of experimental 
studies into the strongly interacting regime, where the optical lattice is strong and the (non- 
rotating) Bose gas is close to the Mott insulating regime, will allow the study of vortex- 
commensurability effects in a regime of strong correlations. 

3.1 Rotating Optical Lattices 

Conceptually, the simplest way to investigate the physics of vortices in a static periodic po- 
tential is to impose a rotating optical lattice on a trapped BEC in an otherwise cylindrically 
symmetric potential well. The optical lattice is at rest in the rotating frame. In the steady 
state, the gas will come to equilibrium in this frame of reference. 

The coupling of a BEC to rotating lattices has been explored in experiments in which 
a rotating mask impose a rotating optical lattice of square or triangular symmetry [134] . 
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Owing to the small angle of the light passing through the mask relative to the rotation 
axis, the periods of these lattices (~ 8/zm) are large compared to the optical wavelength. 
An alternative geometry has been proposed that could allow rotating optical lattices with 
smaller lattice constants, of order half the optical wavelength[26]. 

For a strong optical lattice, in which atoms are restricted to a single Wannier orbital at 
each site |135] 1136] . schemes have been proposed by which to imprint the required Peierls 
phases to simulate the motion of a charged particle in a uniform magnetic field. The uniform 
magnetic field has the same effects as the Coriolis forces of uniform rotation; we shall 
therefore use "rotation" and "magnetic field" interchangeably, assuming that a harmonic 
confinement potential is applied to cancel the centrifugal effects of rotation. One class of 
method makes use of optically induced transitions to other hyperfine levels of the atom 
to imprint geometric potentials on the tunnelling matrix elements |137[ 128) . In another 
proposal, time-dependent potentials are imposed to modulate the lattice potentials and 
tunnelling amplitudes in register. Provided this modulation is sufficiently fast compared 
to the timescales relevant for the physics of the cold atomic gas, these have the effect of 
imprinting the necessary phase modulations to simulate a uniform magnetic field [29J. 

The results of these studies [T3"T1 12"51 12"9"] are proposals for generating systems that are 
well-described by the Bose Hubbard model in a uniform magnetic field 

H = -J^\a\a,e iA ^ +a]aie lA ^ + -U ^ ft^n* - 1) - fi ^ hi (102) 



(id) 

i(t) 



where a\ are the field operators for boson creation/destruction at the lattice site i, and 
hi = d\ai is the number operator. J is the tunnelling energy and U is the on-site boson 
interaction. The ratio size of J/U can be varied by changing the strength of the optical 
lattice[135l 1 136] . The sites i are defined on a 2D lattice, which, for convenience, we shall 
consider to be a Bravais lattice (all sites are equivalent, with one site per unit cell). To 
simulate a uniform magnetic field (or uniform rotation), the link phases Aij — —Aji must 
have the gauge-independent property that their sum around any unit cell of the lattice (the 
discrete version of the line integral § A ■ dr) is a constant 

]T Aij = tor f (103) 

unit— cell 

where < / < 1 is the mean number of vortices (flux quanta) per unit cell of the lattice. 

A related scheme has been proposed in which "non-abelian" gauge fields can be generated 
in optical lattice systems [138] . We will restrict attention to the physics of the above model 
with an abelian gauge field (|102[) . 



3.2 Weakly Interacting Regime: Vortex Pinning 

For a weakly interacting Bose gas, and in the regime where the 2D density of particles is 
sufficiently large, ri2& n v ,n a i tc it is appropriate to treat the system within mean-field 
theory. 

This was first discussed for a uniform system, within an effective model for the vortex 
interactions in the presence of a square lattice pinning potential P3H] • The maxima of the 
potential act as pinning centres for the vortices, owing to the particle depletion at the vortex 
cores. At a vortex density of / = 1 it was found that there exist three phases: the unpinned 
triangular Abrikosov lattice; the strongly pinned square lattice with one vortex per lattice 
site; and an intermediate phase in which half of the vortices are pinned. The phase diagram 
is a function of the strength and period of the pinning potential relative to, respectively, the 
chemical potential and the healing length of the BEC. Extensions to other vortex densities 
/ < 1 and calculations of the collective mode dynamics are reported in Ref. [140] . 
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A study of a rotating square lattice over a wide range of vortex fillings in a harmonically 
trapped atomic BEC is reported in Ref . [141] . A direct numerical solution of the GP equation 
in a trap containing ~ 100 lattice sites was performed. A very rich range of configurations 
was found, including states with doubly quantized vortices pinned to certain lattice sites, 
and the appearance of domain walls between regions of (locally) different order. Owing to 
the many competing phases, the inhomogeneity of a trapped gas leads to a very much more 
complex situation than the uniform case. 

The experiments reported in Ref. [134j studied the effects of rotating triangular [142] and 
square lattices. The locking of a BEC (initially rotating at a different rate to the opti- 
cal lattice) to a rotating triangular lattice was investigated. No clear enhancement in the 
locking was found at the commensurability point / = 1. For a rotating square lattice, the 
structural phase transition was observed from the triangular vortex lattice (at weak pin- 
ning) to a square lattice (strong pinning), at a pinning strength in rough agreement with 
the predictions of Refs.[139] 1141] . Incommensuration effects were found to be difficult to 
observe in experiment [134] . as the finite number of vortices can compress or expand to stay 
close to commensurate conditions. 

Another approach to the physics of this mean-field regime is applicable in the case of 
a very strong optical lattice, for which the Bose-Hubbard model (|102[) can be used. Under 
conditions of weak interactions U <C J and when the mean number of particles per lattice 
site is large. [T43] 1144] one can replace the Hamiltonian (1102j) by the phase representation 

h = - E ^ cos (*< - e i +^)-|Ep ( 104 ) 

where Jjj = I^JniUyJ with m the mean particle number per site. This model is familiar 
from studies of "frustrated Josephson junction arrays" . In the limit U -C J it reduces to 
the classical uniformly frustrated XY model, the groundstates of which have been found for 
infinite uniform lattices at certain rational values of the vortex density / 145, 1461 1147] . The 
groundstate at /= 1/2 is two-fold degenerate. A Gutzwiller mean-field study [3"T] of the 
Bose-Hubbard model in this regime J <C U <C J for / close to 1/2 shows the appearance of 
two competing groundstates. In Ref. |144] the influence on the classical groundstates of the 
inhomogeneity in J due to the harmonic trap has been studied. 

3.3 Strongly Interacting Regime 

Understanding the nature of the groundstate of the Bose-Hubbard model in a magnetic field 
(|102[) poses a very difficult theoretical problem. It involves the interplay of several very 
interesting and complex aspects of physics. 

• As discussed above for weak interactions, there is the interplay of commensurability 
of the vortex density n v with the lattice site density n s i te , controlled by / = n v /n s it e - 
This same commensurability is responsible for a complex form of the energy spectrum 
of a single particle[148j, which has a fractal structure referred to as the "Hofstadter 
butterfly" [159]. 

• In the continuum limit (in the absence of any pinning potential), there can appear 
strongly correlated phases related to FQH states when the 2D particle density U2d is 
comparable to the vortex density n V) controlled by the filling factor v = n2d/n v . (See 

USX!) 

• In the absence of rotation, / = 0, the interplay between interactions U and kinetic 
energy J leads to the physics of the Mott-Hubbard transition for bosons [l35l 1136] . 
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The relevant dimensionless parameters controlling the properties of the system are the 
vortex filling factor /, the particle filling factor n, and the ratio U /J of particle interactions 
to the bandwidth. In addition the properties can depend on the geometry of the lattice on 
which the system is defined. (Unless stated otherwise, it shall be assumed that the lattice 
is square.) 

3.3.1 Mean-Field Theory 

The groundstates of the Hamiltonian (] 1 02[) have been studied within mean-field theory by 
several authors, focusing on different aspects of the very rich physics that can arise. 

The nature of the vortices and vortex lattices close to the Mott-insulating states (n is 
integer) were studied in Ref. |150j . It was found that the vortex cores tend to have a local 
particle density that is equal to that of the nearest Mott state, which can be larger or smaller 
than the average background density. 

In Ref. [143) a mean-field treatment of the quantum phase model (|104[) has been used to 
study the transition at / = 1/2 into the Mott insulating phase. The boundary of the Mott 
phase has been studied in detail in Ref. [151] as a function of the vortex filling /. It was shown 
to follow a dependence that is set by the extremal energy of the corresponding one-particle 
state in the Hofstadter spectrum. That the phase boundary of a strongly interacting phase 
(the Mott insulator) can be described in terms of a single particle property is a surprise; 
the connection arises from the fact that the boundary is determined by the appearance of 
dilute particles, the properties of which are governed by the single-particle spectrum [152| . 

Ref. |152j has studied in detail the vortex lattice structures close to the Mott-insulating 
state. It is shown that, in contrast to the case of weakly interacting Bose gases [1321 H41) 
the vortex cores do not necessarily get pinned to the maxima of the potential (i.e. the 
plaquettes of the lattice) ; vortex lattices are found in which the vortices are pinned to the 
minima of the potential (the sites of the lattice). This is related to the presence of Mott 
phases at vortex cores [HQ], which can act to increase the particle density at the vortex core 
as compared to the average particle density, making it energetically favourable to locate the 
core in the potential minima. 

3.3.2 Strongly Correlated Phases 

Going beyond mean-field theory involves considerations of the effects of quantum fluctuations 
of the vortices. 

The quantum fluctuations of a single vortex have been described in Ref. [153] , for J U, 
when the interaction term in Eqn. (ll04p is small. Fluctuations lead to a finite vortex mass, 
which is predicted to have experimental consequences in Bragg spectroscopy measurements. 

For a dense array of vortices forming a vortex lattice, one can expect that, if sufficiently 
strong, quantum fluctuations will disorder the vortex lattice. A novel form of disordered 
vortex phase was proposed for the "dice" lattice, at vortex density / = l/3 |154j . It is 
argued that in this situation, the first effect of quantum fluctuations of the vortices is to 
lead to a so-called "vortex Peierls" phase: in this phase the vortices are partially delocalized, 
tunnelling between pairs of nearest-neighbour sites, while retaining the broken translational 
symmetry expected for a conventional vortex lattice. 

For very strong quantum fluctuations, one expects the vortex lattice states to disorder 
completely and form strongly correlated liquid states, analogous to the incompressible liquid 
states discussed for rotating Bose gas in the continuum, see ^2.4.31 Indeed, the Bose-Hubbard 
model in a magnetic field (|102p contains all of the physics of the continuum (as discussed in 
U2J in the corner of parameter space with vortex filling and particle filling are small n, / <C 1: 
in this continuum limit, the groundstate depends only on the filling factor v = n/f and 
includes the correlated phases described in ^2.4.31 



42 



The influence of a lattice potential on the bosonic Laughlin state (|T5|) has been studied 
in large-scale exact diagonalization studies on periodic lattice geometries [29l [26] . intended 
to represent the bulk region of a system with a large number of lattice sites. Fixing the 
filling factor (|65|) to v = n/f = 1/2, the influence of the lattice was varied by increasing 
the vortex density / from small values (when a continuum approximation is valid), to large 
values (when the lattice structure is important). It is found that the groundstatc is well 
described by the Laughlin state for sufficiently small /. Calculations of the wavefunction 
overlap indicate that the (continuum) Laughlin wavefunction accurately describes the exact 
groundstate for / < 0.25 29, 155J : a computation of the Chern number associated with the 
topological order show that this takes the value expected for the v — 1/2 Laughlin phase for 
/ ;$ 0.4(261, IT5"5"] . The excitation gap of the Laughlin state is estimated from the numerics as 
a function of / and U/J. In the hard-core limit U ^> J, the gap above the v = 1/2 Laughlin 
state is found to be as large as ~ 0.25J at / ~ 0.11 |26j . 

Exact diagonalization studies have been performed to investigate the properties of small 
clusters of sites with open boundaries, corresponding to rotating lattices in a bounded trap 
geometrv [156[ [157, 158 . As a function of the rotation rate, the groundstate shows transi- 
tions between states of differing "quasi- angular momentum" - the conserved quantum num- 
ber associated with the point-group symmetry of the lattice in the rotating frame. These 
transitions represent the entry of vortices into the system fl 5 7) . 

Numerical constraints for the Bose-Hubbard model limit exact diagonalization studies 
to very small system sizes (typically the number of particles is N < 5 in Rcfs. 29, 26[ I156|, 
11571 I158| ). This makes the investigation of possible strongly correlated phases of bosons 
on rotating lattices very challenging. An analytic proposal for strongly correlated quantum 
Hall states of bosons has been given in Refs.|30[ [31]. The construction is based on the 
observation that the single-particle energy spectrum close to rational fractions / = l/q + 5 
(I, q integer, 5 <g; 1) appears to form a set of q states (which differ strongly on the scale of the 
lattice constant a) , but each of which is modulated by the same Landau level wavefunctions 
over a large lengthscale of order \fl/5 times the lattice constant. For 5 sufficiently small, 
one can treat the system within a continuum theory, but now for q species of bosons. The 
interspecies interactions are determined by the microscopic wavefunctions. At / = 1/2 + 5 
it is argued that in the continuum limit (n, S <C 1) one can construct exact groundstates of 
the model with on-site interaction (|102p ; these are formed from the "(221)" quantum Hall 
state[T5"9"] [see Ean.([ITQ"])] at the effective filling factor v' = n/5 = 2/3. 

The physics of strongly interacting bosons in rotating lattices is a very rich and interesting 
topic, with the potential for exotic strongly correlated phases. Much of the phase diagram 
remains to be explored. It is an area where experiment is likely to raise surprises. 

4 Rotating Multi-Component Bose Gases 

The physics of ultra-cold atomic gases is very much enriched by the possibility to trap and 
cool more than one atomic species. Mixtures of degenerate atomic gases offer a vast array 
of different situations and parameter regimes in which new and interesting phenomena can 
emerge. Consistent with this wide range of parameters, there have been many theoretical 
studies of vortices and vortex lattices in multi-component Bosc gases. Here, I describe only 
those studies that are most closely related to the topics discussed in S}2]for the one-component 
Bose gas. Specifically, I shall focus on results in the 2D LLL regime (|2"TT) appropriate for 
conditions of rapid rotation, and shall limit attention to the two-component Bose gas and 
the spin-1 Bose gas. 
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4.1 Two-Component Bose Gases 

The bulk properties of an ultra-cold two-component Bose gas, with the components labelled 
by i — 1,2, are characterized by the two masses Mj and the set of all mutual two-body 
s-wave scattering lengths, q,-j |160j . A two-component BEC is described by the condensate 
wavef unction 




(105) 



and the mean-field interaction energy is 

E ^ = \ I £ JfcWW d 3 r (106) 

where gtj = 2%Ti a%illHj with fiij = MiMj/(Mi + Mj) the reduced mass |161] . 

For the case of equal interactions gij — g there is an exact symmetry of the interaction 
energy under SU(2) rotations of the condensate wavefunction (|105p . It is then convenient to 
parameterize the condensate wavefunction in a manner that makes this invariance manifest, 
writing 

^ Mr) J ~ V \ cos(0/2)e-^/ 2 

where n is the particle density, \ a global phase, and 9 and 4> parameterize the spinor co- 
ordinates by the polar and azimuthal angles on the surface of the sphere. For SU(2) invariant 
interactions, g^ = g, the interaction energy is invariant under rotations of 6, <fi, as well as 
of the global phase x- F° r experimentally relevant situations, the interaction parameters 
can be close to the SU(2) invariant condition^ Departures from the SU(2) lead to weak 
anisotropies in the energy as a function of and oi> )165| . 

For a uniform 3D system, stability of the individual condensates requires 311,(722 > 0. 
Within this stable regime, depending on the interactions the two components can be miscible 
or immiscible. The condition for immiscibility, in which case the BEC phase separates into 
spatially separated regions of component- 1 and component-2, is [1601 1166] 

3ii522 < 9I2 ■ (108) 

Two-component BECs have been formed from the trapping of two hyperfine states of the 
same atomic species (for 87 Rb |167] . and for 23 Na [168] ). and from the trapping of mixtures of 
two different atomic species (for 41 K- 87 Rb[169l H2Q] and 85 Rb- 87 Rb[17T])- Experiments on 
these two-component condensates show evidence for both miscible and immiscible regimes 
compatible with (|108p . and have investigated the transition between these regimes as the 
scattering lengths are varied fTTOl 1171] . 

When the two components are two hyperfine states of a single atomic species, the masses 
are identical Mi = M 2 = M; when they correspond to two different atomic species, then 
Mj 7^ M 2 . The rotational properties of the two-component gas are very different in these 
two cases, which we therefore discuss in turn. 

4.1.1 Rotating two-component BEC: Equal masses 

Single Vortex States 

Experimental studies of quantized vortices in atomic BECs were initiated by experiments 
on a two-component BEC |172] , Coherent excitation of a 87 Rb BEC in a single hyperfine 

12 For the \F = l,mp = —1) and \F = 2,mp = 1) states of 87 Rb, the ratios are gn : 922 '■ 912 = 
1.024 : 0.973 : 1.0 [T62l IT63"] . For the \F = l,m F = 0) and \F = l,m F = 1} states of 23 Na these are 
911 : 922 : S12 = 1 : 1-035 : 1.035HM]- 



(107) 
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state (which could be either one of the two trapped states) was performed in such a way 
as to convert atoms into the other trapped hyperfine level with one unit of circulation. 
The imprinted configuration therefore had a single quantized vortex for the particles in 
the excited component. The phase of the rotating component was directly imaged in an 
interference technique, to confirm the quantized circulation. The subsequent evolution of 
the two-component system was studied, showing evidence for differences in the dynamics 
and stability depending on which of the two hyperfine levels was excited into a circulating 
state. 

The structure of the axisymmetric vortex state created in these experiments is easily 
written in the 2D LLL regime, for which the condensate wavefunction takes the form [173] 

( £ ) - ( k ) * ■ <™> 

where a and b are amplitudes whose ratio determines the relative populations of the two 
components iVi/A^. This wavefunction describes a state in which component-2 has one 
unit of angular momentum while component- 1 is at rest. Note that, despite the density 
suppression in the core of the vortex of component-2, there remains non-zero particle density 
for component-1 in the core region. This is a general feature of the vortices in two-component 
BECs with repulsive interspecies interactions g 12 > |160U174l 1175111761 117311165] . Since the 
density does not vanish, one can use the representation (|107[) to describe this vortex as a spin- 
texture in the variables (0, 4>, x). In terms of this description, the vortex configuration (| 1 09|) 
has a non-trivial topological classification. It is referred to as a "half quantum vortex" |176j 
or (two-dimensional) "Skyrmion" |173j . The properties of these spin textures have been 
investigated in detail for interactions that (weakly) break the SU(2) symmetry[165j. The 
results show the appearance also of vortex states which break the axial symmetry of the 
trap. 

Beyond mean-field theory, there exist exact results for the case of SU(2) invariant 
interactions, — g, in the 2D LLL limit. For small values of angular momentum, 
L < min(iVi, N2), the results of exact diagonalization studies on N = N\ + N2 < 8 particles 
have been used to motivate a conjecture for the exact (analytic) groundstates [T77] . The nu- 
merical studies show that the groundstate contains only single particle states with angular 
momenta m < 1; such states can be shown to be exact eigenstates of the SU(2) interac- 
tions. Similar conclusions hold also for the rotating spin-1 Bose gas with SU(3) invariant 
interactions [1 78] . discussed in £14.21 

Many Vortex States 

Theoretical studies of lattices of many vortices were performed by Mueller and Ho |179j 
within the mean-field LLL regimeF^I As for the case of the one-component BEC, in the LLL 
limit the condensate wavefunction may be expressed in terms of the positions of the vortices. 
However, for the two-component condensate wavefunction there is a set of vortex locations 
for each of the two components. For the case gn = 322 = <? 7^ 9vz studied in Ref. |179j the 
two components are identical, so one expects each to have the same vortex structure up to a 
relative translation. The groundstates were determined as a function of 512 / g for an infinite 
lattice. A series of vortex lattice phases was found, with structures depicted in Fig. [TBI 
These results illustrate the very rich physics that is possible in a rotating two-component 
BEC. For gyijg > 1, corresponding to the regime of phase separation (|108| for the non- 
rotating gas, the groundstate consists of interlaced rectangular arrays of vortices. For large 
gnlg the vortices become closely spaced along one direction, leading to vortex sheets; the 
particle density then tends towards an array of parallel stripes in which one component has 
low density and the other a high density. 

The vortex lattice states have been studied beyond the LLL and including the effects of 
the trap |180l 1181] . The results show qualitatively similar phases to those described above. 

13 Within mean-field theory the results are identical for the quasi-2D and 3D LLL regimes. 
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Figure 15: Vortex lattices in the mean- field LLL regime of a two-component BEC for 
.9n = .922 = .9, and at different values of gvijg. The filled and shaded circles show the 
locations of vortices in each of the components. For .912 /<? < the interspecies interactions 
are attractive, so the vortex cores of the two lattices attract each other. The vortex lattices 
of the two components each form a triangular vortex lattice, with coincident vortices (a). For 
912/9 > the interspecies interactions are repulsive, so the vortex cores repel each other. 
A series of lattices is found which appear as interlaced triangular lattices (b); interlaced 
rhombic lattices (c); and finally interlaced square (d) or rectangular (e) lattices. At the 
SU(2) point (712/3 = 1, the interlaced rectangular lattices of the two components can be 
viewed as a triangular lattice if the differences of the two components is ignored. The final 
panel illustrates the variational parameters used in the study. [Reprinted figure from: E.J. 
Mueller and T.L. Ho, Phys. Rev. Lett 88, 180403 (2002). Copyright (2002) by the American 
Physical Society. ] 
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Stripe- like states (vortex sheets) appear in the regime of phase separation. Dynamical modes 
have been studied in the uniform system in the LLL [182] and beyond the LLL in the trap 
geometry [HI]. 

Experimental studies of the vortex lattices in rotating two-component 87 Rb BECs have 
been reported in Ref. [183] . The interspecies interactions are repulsive gi2 > so the vortices 
in the two components repel each other. The interactions are close to the SU(2) invariant 
values gn ~ gi2 — 322 (see footnote 1 on page [44]), but just inside the regime of phase 
separation (]108|) . The experimental results showed evidence for the appearance of a stable 
vortex lattice phase with two inter-penetrating lattices of approximate square symmetry, 
consistent with the phases found in Refs. |179[ I180j . 

Exact results for the interacting two-component system with many vortices are very 
limited. 

In the 2D LLL regime, the groundstate at high angular momentum can be determined 
analytically for repulsive contact interactions. Recall that, for spinless bosons with repulsive 
contact interactions, the exact groundstate at high angular momentum is the Laughlin state 
([79]) with total angular momentum L — N(N — 1) corresponding to filling factor v = 1/2. 
For the two-component case with g^ > 0, and for equal numbers of atoms in the two 
components N\ = N2 = N/2, the exact zero-energy groundstate is an incompressible liquid 
described by the "(221)" bosonic Halperin state |159] 

i<j i<j i,j 

where ^ are the co-ordinates of the N/2 particles in component-1 and ^ the co-ordinates 
of the N/2 particles in component-2. This state can be viewed as a fluid of pairs of particles 
forming spin singlets. This function vanishes when the co-ordinates of any two particles 
coincide, so it is a zero-energy eigenstate of the contact interaction Hamiltonian; it is the 
groundstate provided the interactions are repulsive g%j > 0. The total angular momentum 
is L = 3N 2 /4 — 1; which corresponds to a uniform fluid state with a total filling factor 
v\ + V2 = 2/3, therefore denser than the Laughlin state. 

A series of non-abelian states that can describe two-component Bose systems in the 
2D LLL was proposed in Ref. [184 j . These "non-abelian spin-singlet" (NASS) states form 
a sequence labelled by k = 1,2,..., and are the exact groundstates of the two-component 
Bose gas for a k + 1-body contact interaction ([55]) . at filling factors v — 2fc/30 For k = 1 
the state in this sequence is the above Halperin (221) state; for k > 2 the NASS states 
are non-abelian phases which are spin-singlet generalizations of the Moore- Read and Read- 
Rezayi states. Owing to the importance of the Moore-Read and Read-Rezayi phases in 
describing the phases of rotating spinless bosons with realistic two-body interactions (see 
§2.4| , it would be interesting to explore the relevance of NASS states for a two-component 
Bose gas with realistic two-body interactions. 



4.1.2 Rotating two-component BEC: Unequal masses 

The rotational properties of a two-component BEC formed from two atomic species raises 
additional interesting features. Since the atomic masses differ, Mi ^ M 2 , if the two compo- 
nents were to rotate (on average) at the same rotation frequency f2, then by the Feynman 
criterion (|10[) the mean vortex densities in the two components would be unequal 

^ = ^ (111) 



'The more general form is v = 2k/(2kq + 3) with q an even/odd integer for bosons / fermions |184| . 
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This imbalance of the average vortex densities competes with an energetic gain from inter- 
species interactions which is obtained by making the vortex densities equal (or commensu- 
rate). 

The multi- vortex states of a two-component BEC with unequal masses have been studied 
theoretically in Ref. 185J. The interspecies interaction is assumed to be attractive gn < 0, 
such that the cores of the vortices in the two components attract. 

For strong attractive interactions, the vortices pair, and the system forms a locked state 
in which the average vortex densities are equal. How does one reconcile the equality of the 
vortex densities in this phase with the above relation (|llip ? Consider a system that has 
come to equilibrium at a drive frequency f2, and that forms a locked vortex lattice of pairs 
of vortices. The locations of the vortex cores must rotate uniformly at the drive frequency 
f2, and they will have some density n v that is determined by the energetics (including 
the inter- vortex interactions). Since the vortex density is fixed to n v , the individual fluid 
components are flowing with average flow fields that lead to rotation of the two fluids with 
rates = hn v / (2Mi) . Clearly Hi ^ and, in addition, neither f^i nor need match 
the drive frequency Q. Thus, there is a relative motion of the fluids with respect to the 
vortex cores. This relative motion leads to Magnus forces on the vortices. One reconciles 
these features by noting that the Magnus forces are exactly balanced by the inter-particle 
interactions that act to pin the vortices to each other [185] . as they must be for consistency 
with the statement that the vortex lattice is at equilibrium at the drive frequency O. 

For weak attractive interspecies interactions the vortices separate and the vortex densities 
are set by those expected for independent condensates, with nip ^ rvp The transition 

between these two regimes is controlled by the parameter \gi2\n/(hfl) where n is the mean 
particle density. Thus, in a trap geometry, the inhomogeneity of density leads to spatial 
variation: a region of locked vortex lattice at centre of the trap, and an unlocked region at 
the outer edge |185j . 

Related issues arise in models of rotating Bose gases close to a Feshbach resonance [TT^]. 
The atoms and molecules (formed from resonant association of pairs of atoms) have unequal 
masses M and 2M, leading to different mean vortex densities for the molecular and atomic 
condensates [186]. It is found that the coherent coupling can lead to the formation of states 
in which the atomic and molecular vortices form units resembling a carbon dioxide molecule. 
The arrangement of the orientations of these units within an approximately triangular array 
leads to many near degenerate configurations and glass-like behaviour. A range of other 
interesting vortex lattice phases is identified [186]. 



Spinor atomic gases may be realized by trapping a high spin atom in an optical trap, 
which does not lift the spin degeneracy. The inter-particle interactions are spin-rotationally 
invariant, and for spin S = 1 must take the form |187l 1188] 



where cq and ci are the only two interaction parameters. In terms of the two s-wave 
scattering lengths, denoted a® and al for scattering in the two scattering channels with total 
spin S — 0, 2, the coupling constants in Eqn. (|112|) are cq — (go + 2<?2)/3, C2 = (32 — 5o)/3, 
where g s = 4nh 2 af /M§ 

For C2 < 0, the system acts to maximize the total spin, which is achieved by aligning 
all spins; this is referred to as the "ferromagnetic" state. It is characterized by a non-zero 

15 Since the components are different hyperfine states of the same atomic species, the mass M is the same 
for all components. 



4.2 Spin-1 Bose Gas 




(112) 
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value of the expectation value of the local spin density (S). The breaking of spin-rotational 
symmetry leads to a manifold of groundstates related by SO(3) transformations|188]. For 
c 2 > 0, the system acts to minimize the total spin; this is the "polar" state. It has vanishing 
spin density, (S) = 0. Still, it breaks spin-rotational invariance, with a groundstate manifold 
described by S 1 x S 2 / Z-? ^18S\ 1189] . Thus, in contrast to the two-component condensate, 
which is characterized by a global phase and a local spinor direction (|107[) . the topology of 
the spin-1 Bose gas is much richer. 

Both ferromagnetic and polar spin-1 BECs have been studied experimentally, with 87 Rb 
(c 2 < 0)[1"90] and with 23 Na (c 2 > 0) [T9T 1 fTM ] [168] . In both of these experimental systems, 
the spin-dependent interaction is weak, with | C2 / Co | *C 1. Consequently, in spin-textured 
configurations such as vortex lattices, it can be energetically favourable for the condensate 
to have mixed character, with regions in which the local order parameter is of ferromagnetic 
character (with non-zero average spin) and regions where it is of polar character (with 
vanishing average spin). 

Single Vortex States 

The resulting textures for single vortices of a spin-1 BEC have been studied theoretically 
within mean- field theory [19^1 11931 [194] , Within the LLL mean-field regime the description 
of these states simplifies [1781 [173]. However, a wide range of textures can still appear, with 
forms that depend on c 2 /co and on the total angular momentum. The angular momentum 
is carried by spin-textures, which can be of the form of "skyrmions" of the ferromagnetic 
order or of "7r-disclinations" of the polar phase. 

In the 2D LLL regime, the exact groundstates at small angular momentum, L < N, have 
been found for the case c 2 = in which case the interactions have a full SU(3) symmetry [178] , 
There exist exact eigenstates of the Hamiltonian for which the only occupied single particle 
states are those with angular momenta m < 1 (for < L < 2N/3), or those with m < 2 
(for 27V/3 < L < N). Exact diagonalization studies confirm the states to be the exact 
groundstates for all system sizes studied. 

Many Vortex States 

The vortex lattices of a uniform S = 1 Bose gas have been studied within mean-field 
theory for the polar regime c 2 > 0[195j. A rich series of phases of five different lattices of tex- 
tured vortices was found, depending on c 2 /co- The vortex lattices were studied within LLL 
mean- field theory in Ref. [178j for both positive and negative c 2 . Eight different symmetries 
of vortex lattice were found, which are illustrated in Fig. 1161 The phases involve arrange- 
ments of various spin-textures of the spinor order, with lattices of Skyrmions, 7r-disclinations 
and conventional vortices [178| . 

Turning to consider rapidly rotating gases beyond mean-field theory, there exist some in- 
teresting exact results. In the 2D LLL regime appropriate for rapid rotation, the groundstate 
at high rotation rate can again be determined analytically for repulsive contact interactions, 
50j<?2 > 0. It is an incompressible liquid state that may be viewed as a generalization of 
the Halperin (221) state, Eqn. (|110p . The exact (zero-energy) groundstate can be viewed as 
a liquid of clusters of triplets of particles in the three spin projections of spin-1 [196, 197 , 
and is the SU(4)|i state in the sequence of Ref. |198j . For N divisible by three (such that 
all particles can form complete clusters) the angular momentum is L = N(2N/3 — 1), which 
corresponds to a filling factor v — 3/4. These clustered states may be generalized to an 
arbitrary number of components n; they form the exact groundstates of repulsive contact 
interactions at filling factors v = n/(n + 1) |197[ [T98] , 

A series of generalized Read-Rezayi states has been proposed for rotating S — 1 bosons [198]. 
The states in this sequence, labelled SU(4)|fc, have the property that their wavefunctions 
vanish when the co-ordinates of k + 1 particles (in any spin component) coincide; they are 
therefore exact groundstates of repulsive k + 1-body contact interactions (|85[) between the 
atoms. They describe incompressible liquid phases at v = involving clusters of 3k 

16 The more general forms are v = 3k/(3kq + 4) for SU(4)|j. and v = k/(kq + 1) for SO(5)|fc with q an 
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Figure 16: (Colour.) Vortex lattices in the mean-field LLL regime for a spin-1 BEC at 
different values of ca/ 'co- The first panel shows the Hammer- Aitoff projection of the colours 
on the spin-sphere, in which the intensity encodes the size of the spin \{S)\. The images 
are the spatial distributions of the expectation values of the local spin density for the vor- 
tex lattices at c 2 /c = -0.1,-0.05,0.01,0.016,0.04,0.1,0.54,0.7,0.9. In the last picture, 
C2/C0 = 0.9, the spin vanishes, and so the intensity shows the particle density, which forms a 
triangular Abrikosov lattice. [Reprinted figure from: J.W. Reijnders, F.J.M. van Lankvelt, 
K. Schoutens, and N. Read, Phys. Rev. A 69, 023612 (2004). Copyright (2004) by the 
American Physical Society. ] 

particles. For k = 1 this is the clustered state |196] I197| I198j described above, which is the 
exact groundstate for two-body contact repulsion between all particles. For k > 1 the phases 
have quasiparticle excitations which obey non-abelian exchange statistics. Although not the 
exact wavefunctions for two-body interactions, these represent interesting trial wavefunctions 
for possible quantum Hall phases of rotating spin-1 Bose gas at v > 3/4. Refs. |198l I178j 
report the results of numerical studies of the exact groundstates in this regime. A large 
overlap of the exact groundstate (of the two-body interactions) with the SU(4)|2 state is 
reported at v = 3/2, which is evidence of the appearance of this non-abelian phase. 

A second generalization of the Read- Rezayi states, labelled SO(5)|fe, and describing S = 1 
bosons at filling factor v = k is proposed in Ref.[198 . These states are exact zero-energy 
eigenstates of a k + 1-body contact interaction, for a special choice of c^Icq. They involve 
clustering of groups of 2k particles. These phases have the additional interesting feature of 
involving fractionalization both of atom number or "charge" (as is common in FQH states) 
and also of spin, their quasiparticle excitations having atom number 1/2 and spin 1/2. 

5 Rapidly Rotating Fermi Gases 

The possibility to cool Fermi gases to quantum degeneracy T < Tf |199] and to achieve 
regimes of strongly correlated quantum phases of two-component Fermi gases[4], raises in- 
teresting questions concerning the properties of degenerate Fermi gases under rotation. 

For a one-component Fermi gas s-wave interactions vanish due to the Pauli principle. 
For an ultra-cold atomic gas (with low relative momentum tik between colliding particles) 
the dominant interactions are in the p-wave channel. The cross-section is typically very 
small, vanishing as <7\ ~ V 2 k 4 as k — ► where the "scattering volume" V is set by the inter- 
atomic forces and therefore typically V <C 1/fc 3 for a cold gas. Nevertheless, it is possible to 

even/odd integer for bosons/fermions 198 . 
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introduce significant interactions between the identical fermions. This can be achieved by 
using atoms or molecules that have strong dipolar interactions|58j : or by use of a "Feshbach 
resonance" to enhance the p-wave interactions [200]. 

For a two-component Fermi gas, there can be s-wave interactions between the two compo- 
nents. Moreover, by use of a Feshbach resonance, these interactions can be made very large. 
In this way, experiments have been able to study the evolution between the BCS paired 
state (for weak attraction between the two species) and the BEC state (of small molecules 
formed from tightly bound pairs of the atoms), through a regime of strong correlations [4]. 

In this section we discuss theoretical predictions for the effects of rapid rotation on one- 
and two-component Fermi gases. 



5.0.1 Non-Interacting Fermi Gas 

The single particle states for a fcrmion in an axisymmetric harmonic trap are the same as 
those discussed in ij2.ll However, we shall now consider states beyond the lowest Landau 
level, so need to generalize the discussion. 

Working in the rotating frame, at angular frequency Q, the single particle energies of 
particle in an axisymmetric harmonic well are 

£ S J .,m,n„ =Ku ± {2n ± + \m\ + l)-nQ.m + nu\\{n\\ +1/2) (113) 

where nj_ = 0, 1, 2 . . . and to — 0, ±1, ±2 are the radial and angular momentum quantum 
numbers for motion in the plane perpendicular to the rotation axis, and riy = 0, 1, 2 . . . is 
the subband index for motion parallel to the rotation axis. If there are additional internal 
degrees of freedom (hyperfine levels), there is one such orbital state for each level. 

In general, one should take account of confinement in all directions, owing to the inho- 
mogeneity of the trapped gas. However, as in the case of bosons, it is useful to consider the 
limit £1 ~ lu±_ in which the system becomes locally homogeneous in the plane perpendicular 
to the rotation axis. The spectrum can then be written 

£^ m>n|| = 2hoj ± (n + 1) + ^|| (n,| + 1/2) (114) 

where now n = 0, 1, ... is the Landau level index and the angular momentum takes the values 
m = — n, —n + 1, — n + 2. [The relationship is obtained by noting n = n± + Qm\ — m)/2.] 
Note that the states in a single Landau level (labelled by to) are degenerate. The particle 
density sets the highest occupied Landau level (maximum n) and subband (maximum ny). 

The properties of the gas depend on the aspect ratio of the trap, u>±/u>i\ . It is of interest 
to consider the two limiting cases: the quasi-2D regime, u>\\ ^> and the 3D regime, 

In the quasi-2D regime, the particle density is sufficiently small or the axial confinement 
is sufficiently strong that only n\\ = is occupied, i.e. /i < Eqq l . One can then define the 
filling factor as before, and as is conventional in the FQHE systems, by 

n 2 d n 2 dh , . 

v = = . 115 

n v 2Mw± y ' 

In this quasi-2D regime, even a non-interacting Fermi gas shows a series of incompressible 
states. These arise due to the fact that the spectrum is discrete, with a set of highly 
degenerate states (over to) at E^^ . These appear each time v is an integer and a Landau 
level is filled. (At this point the chemical potential has a discontinuity equal to the Landau 
level spacing, 2ftw_L.) If there are additional internal degrees of freedom these incompressible 
states occur at those filling factors for which atoms of all species fill an integer number of 
Landau levels. 
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The 3D regime ([33)) is relevant for anisotropic traps uu <C uj_l- Here, for each Landau 
level state, the motion along the rotation axis can be viewed as a quasi-lD gas, and the 
energy is more conveniently written as 

h 2 k 2 

= 2f l w ± (n+l) + ^l ) (116) 

where fen is the wavevector along the rotation axis. In this case, the system does not show 
incompressibility. However, there is a \j\fE divergence in the density of states at the 
position of each Landau level, associated with the quasi ID motion. 

The existence of Landau level quantization is predicted to lead to signatures in the 
density distribution for a trapped Fermi gas, when the rotation frequency is close to the 
trap frequency fl ~ u> i [201] . The experimental signatures depend on the aspect ratio of 
the trap: they are strongest for the quasi-2D limit, where there is strict incompressibility 
leading to plateaus in the density distribution; but features survive also in the 3D case due 
to the singularities in the density of states [2171] . 



5.1 One Component Fermi Gas 

Under conditions of rapid rotation, and at sufficiently low particle density, the Fermi gas 
will be in the quasi-2D LLL regime (n — n\\ =0). At the single particle level, the system is 
then degenerate when the filling factor is non- integer, v < 1. The groundstate is determined 
by minimizing the inter-particle interactions. 

In the following I describe the results of studies of the effects of two important forms of 
inter-particle interactions on the properties of a rapidly rotating one component Fermi gas. 



5.1.1 Dipolar Interactions 

The properties of a Fermi gas in the 2D LLL have been investigated for the case of dipolar 
interactions. Antisymmetry of the wavefunction requires the relative angular momentum of 
any two particles to be odd, so the inter-particle interaction is characterized by the Haldane 
pseudo-potentials ([5Tj) at odd m. The relevant energy scale can be large for fermionic 
molecules with electric dipole interactions 202 , leading to the possibility of stable correlated 
phases in rapidly rotating dipolar Fermi gases. 

The possibility to use dipolar interactions to stabilize the v — 1/3 Laughlin state 
was described in Ref. [202] . The energy gap for creation of a quasi- hole excitation (which 
contributes to the incompressibility of the state, see ^2.4.3[) was estimated to be A q h = 
(0.9271 ± 0M9)C d /£ 3 for a trap with spherical symmetry {w x = ^||)[2Q2]. 

The groundstates of small numbers of dipolar-interacting fermions (N < 10) in the 
2D LLL in a harmonic trap have been investigated numerically[203j over a wide range of 
different rotation rates (angular momenta). The results were analysed within the framework 
of the composite fermion construction, in which the interacting Fermi system is represented 
by non- interacting composite fermions [841 [21]. This construction works well for analogous 
studies of harmonically trapped electrons in the FQHE regime (i.e. Coulomb forces instead 
of dipolar forces) [86] . As described in §2.3.31 a generalization of the composite fermion 
method to treat rotating bosons in harmonic traps also works well[83]. For the case of 
dipolar interacting fermions, some discrepancies are reported from the apparent predictions 
of CF theory [503]. However, these appear to arise from an incomplete application of the 
composite fermion theory for some of the groundstates found in numerics To check the 

17 Most of the discrepancies reported in Ref. 203 arise from the neglect of a set of the "compact" CF states 
that account for the physics in other situations 86 83 . Specifically, for N = 10, the CF theory predicts a 
set of compact states at L = 73, 80, 85, 93, 95; these values fill out most of the missing entries in Table (4) 
of Ref. _203_ . The results reported therefore appear consistent with CF theory describing the high angular 
momentum states accurately. 
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accuracy of the CF description of these states, or to determine whether the states are of a 
different character [203J , it would be useful to compute overlaps of the exact groundstates 
with the trial composite fermion states [86| [83]. 

The fact that dipolar interactions are long-ranged allows the stabilization of a Wigner 
crystal state at very low filling factor v <C 1. The transition from Laughlin states to the 
Wigner crystal phase is predicted to occur at v < 1/7[123 , close to the value at which the 
transition is believed to occur in FQHE systems. 

In the case of a quasi-2D fermionic gas with filling factor larger than 1, v > 1 — i.e. when 
more than one Landau level is occupied - the nature of the groundstate of the particles in 
the partially occupied Landau level (let us label this level by n) are determined solely by the 
interactions. However, now the Haldane pseudo-potentials should be computed in this nth 
Landau level. Consequently, the groundstate at a filling factor v = n + v' (0 < v' < 1) can 
differ from that at v = v' . Studies in the quantum Hall literature show that Hartree-Fock 
mean- field theory is accurate when n 3> 1 [65] . For large filling factors, the mean-field states 
of dipolar atomic Fermi gases are expected[67] to be qualitatively similar to those of electrons 
in semiconductor systems [64 | 1651 . Depending on the filling of the partially occupied Landau 
level, the groundstate will be a "stripe" or "bubble" phase[64, 65j of similar translational 
symmetry to the mean-field states discussed for bosons in ^2.2.31 

5.1.2 p-wave Feshbach resonance 

Although weak at low collision energy, identical fcrmions experience non-zero p-wave in- 
teractions. These interactions can be resonantly enhanced in the vicinity of a Feshbach 
resonance - i.e. close to the threshold for formation of a new bound state in the p-wave 
channel. Such resonances have been studied in ultra-cold gases of 40 K[200, 2041 1205] and 
6 Li [2TMl2irni2l?g] . 

Close to a p-wave resonance, at low collision energy [when the wavevector max(/ci?, l/£) 
is sufficiently small] the p-wave scattering phase shift takes the form 

6i(k) = (l/3)fc 3 a?. (117) 

This defines the scattering length a\. The quantity a\ is often referred to as the "scattering 
volume" [209J , the modulus of which diverges at the position of the resonance. 

In the lowest Landau level, the interaction may be represented by the pseudo-potentials. 
For the p-wave interaction, the only non-zero pseudo-potential is in the m = 1 channel, 
with[210] 

[2h 2 a 3 , s 

V > = ^M^- < 118 > 

This two-body interaction is the special case for which the v = 1/3 Laughlin state is 
exact [2"TT] . For v < 1/3 the groundstate is highly degenerate, and may be viewed as a 
gas of 1/3 anyons of the Laughlin state. 

Analytic results do not exist for v > 1/30 or for the excitation spectrum above 
the v = 1/3 groundstate. Numerical studies on the sphere [2T0] have investigated the 
groundstates and low-energy excitations above the Laughlin v = 1/3 state and the "Jain 
sequence" [84] states at v = 2/5,3/7. These studies show that in the thermodynamic limit 
N — ► oo the excitation gaps appear to converge to non-zero values, establishing these states 
as incompressible and determining numerical estimates values for the excitation gaps. The 
overall energy scale is set by the V\ pseudo-potential. Experimental observation of these 
states will require cooling to temperatures less than ~ Vi/ks- The issue of stability of 
a degenerate one-component Fermi gas close to a p-wave resonance |2051 12121 [208 is an 
interesting open question, as is the effect of rapid rotation on this stability. 

18 One exception is the trivial statement that v = 1 is a fully filled Landau level. 
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Figure 17: Schematic diagram of the variation of the s-wave scattering length close to a 
Feshbach resonance as a function of the detuning AB = (B — B rcs ). 

5.2 Two-Component Fermi gas 
5.2.1 BEC-BCS crossover 

Degenerate gases of two species of fermions may be formed either from two atomic species, 
or from two hyperfine states of the same fermionic atom [2l51 1214) . I shall denote the two 
components as spin-up and spin-down. For simplicity I shall assume that the two components 
have equal densities and equal masses M. 

Since intra-species s-wave interactions are excluded by the Pauli principle, the dominant 
interactions at low energy are the inter-species s-wave interactions. By use of a Feshbach 
resonance, these interactions can be made strong. A Feshbach resonance arises when, close 
to the energy of the two dissociated atoms, there exists a bound molecular state of the 
two atoms in some different hyperfine levels. Owing to the different hyperfine structure of 
the bound state from the dissociated atoms, by varying a magnetic field one can cause the 
energy of the bound state to cross the energy of the two dissociated atoms. Hybridization 
of the molecular level with the unbound atoms leads to a Feshbach resonance in the s-wave 
scattering amplitude. 

For a detailed discussion of the physics of the resonant scattering, the reader is referred 
to Ref. [4] . In brief, the consequence of the formation of the bound state is that the s-wave 
scattering length a s passes through a divergence, 



as illustrated in Fig. 1171 On approaching the resonance from the side on which the bound 
state has not formed (B > B ICS ), the scattering length becomes large and negative; ap- 
proaching the resonance from the side on which the bound state has formed (B < B res ), the 
scattering length becomes large and positive. For short-range interactions, the properties 
of the (balanced) two-component Fermi gas close to the Feshbach resonance depend on the 
parameter kpa s . [We denote the Fermi energy ep = h 2 k F /(2AI).] 

The nature of the groundstate can be understood in the limiting cases far from the 
resonance, when fci?|a s | is small. 

When a s is small and positive, there exists a bound-state even for one pair of atoms. 
The groundstate of a many particle system consists of tightly bound pairs of atoms of the 
two species; the resulting small bosonic molecules behave as a weakly interacting Bose gas 
and form a Bose-Einstein condensate at low temperatures. This is the "BEC side" of the 
resonance. 

When a s is small and negative, the weak attractive interaction causes a BCS instability 
of the Fermi surface, leading to a BCS paired superfluid. There is an energy gap that, in 



a s (B) oc 



1 
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this BCS limit of weak interaction, is exponentially small 

a ~~^h*7Ki)- < 120 » 

This is the "BCS side" of the resonance. 

As the resonance is swept, such that l/kpa s varies from the BCS to BEC regime, it 
is believed that the groundstate of the Fermi gas evolves smoothly. This "BEC to BCS 
crossover" has been confirmed in recent experimental studies of several groups [3]. 

Cold atomic gases allow additional novel physics to be studied by imposing an imbalance 
between the populations of the two components. In conventional superconducting systems, 
a population imbalance can be achieved by using the Zeeman splitting in a magnetic field. 
However, a magnetic field also induces orbital effects, precluding a study of the effects of 
the density imbalance alone. In atomic gases the density imbalance can be achieved without 
any orbital effects [2T5I I2T5] . 

New physics can also emerge in atomic gases in the complementary situation: that is, 
to study purely the orbital effects of a magnetic field, without inducing any population 
imbalance. This can be achieved by setting the densities of the two hyperfine levels to be 
equal, and introducing the orbital effects of the magnetic field by rotation. This is the 
situation that shall be described in this section. 

The effects of moderate rotation are well understood. At all points across the crossover, 
the low energy phase is a condensate, and responds as a superfluid of pairs of fermions. 
Rotation causes the formation of an array of vortices with density, 2(2M)f2//i, which is the 
Feynman result (|10p with the mass, 2M, of the pair of particles. The formation of a lattice of 
quantized vortices has been demonstrated in experiments described in Ref. [217j . The results 
reveal the superfluid nature of the system across the whole range of the crossover regime, 
and are consistent with vortices having a quantum of circulation of h/(2M) associated with 
condensation of pairs of fermions. 



5.2.2 Effects of Rapid Rotation 

One expects the effects of rapid rotation to be rather different in the BEC and BCS sides 
of the resonance. 

Far on the BEC side, the effects of rapid rotation are similar to those discussed above for 
rapidly rotating bosons. In the quasi-2D regime, the mean-field vortex lattice will survive 
up to vortex densities at which the filling factor ([57)) is of order one, at which point there is 
a transition to strongly correlated states. In the 3D regime, the transition to the strongly 
correlated regime occurs at high vortex density according to Eqn. ([M|) . 

Far on the BCS side, the effects of rotation are similar to the orbital effects of magnetic 
field on type-II superconductors. For conventional superconductors, the upper critical field 
is well described by the semi-classical approximation of Gorkov|218, 219, 220J. Translating 
to the quantities of interest to us, this predicts the transition to a normal state wherT^I 

hn> 4.24A 2 /ef- (121) 

One can understand this semi-classical result in terms of a balance of energies. On the 
one hand the condensation energy per particle on forming the superfluid is of order A 2 /cf- 
This should be compared with the kinetic energy of the superfluid flow. Since the su- 
perfluid cannot rotate as a rigid body, in the rotating frame of reference the superfluid 
has a non-zero kinetic energy. To estimate this kinetic energy, note that the vortex spac- 
ing a v ~ ^ 2(2M)o j-1/2 sets a typj ca j velocity h/Ma v , so the kinetic energy per pair is 

19 Note that, in the BCS limit, this is a much lower rotation rate than that required to put one flux 
quantum through one Cooper pair, which would be hfl ;J A. 
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M{h/Ma v ) 2 ~ h^l. The semi-classical formula states that once this kinetic energy is larger 
than the condensation energy, it does not pay to form the superfluid, and a rigidly rotating 
Fermi gas has lower energy. 

How the transition between the BEC and BCS regimes occurs in a rapidly rotating gas 
is an interesting question. 

The mean-field theory for a rotating Fermi gas in a harmonic trap has been studied 
in Ref. [221j across the whole crossover regime. This theory gives results on the BCS and 
BEC sides that match with the above expectations in these limits, and provides a useful 
way to connect between these two regimes. The approach works explicitly with a trapped 
gas, and the mean-field theory has the benefit of being accurate in the limit of a "narrow" 
Feshbach resonance [222]. The narrow resonance limit has the feature that the molecular 
bosons that form have vanishing interactions. As a result, for a trapped gas at any Q =/= luj_ 
the molecular gas on the BEC side is confined to the lowest harmonic oscillator state and 
carries zero angular momentum. This is an unphysical situation for rotating atomic Fermi 
gases which interact, typically, through "wide" Feshbach resonances. However, it is unclear 
under what conditions (if any) this might affect the validity of the predictions for rotating 
Fermi gases beyond the narrow resonance limit. 

BCS Regime: Upper Critical Rotation Rate 

A complementary approach has been described in Ref. [2231 1224] . These works study 
mean-field theory in a situation in which the trapping potential is weak, so the system is in 
the uniform 3D regime, 

fl ~ w±_ > cj|| , (122) 

with single particle energy spectrum (|116|) . The linearized BCS mean-field equations were 
used to determine the critical temperature below which the system is unstable to the for- 
mation of a BCS superfluid, as a function of the rotation rate tt. 

In Ref. [223 . , the effects of a trap are introduced within the local density approximation 
(see §2.5]) . This gives rise to inhomogeneity in the trap, with normal and superfluid regions 
co-existing at different radii. An upper critical rotation rate is determined as a function of 
detuning, at which the entire cloud becomes normal for a temperature of ksT/hfl = 10~ 3 . 

In Ref. [224] . similar mean-field equations were studied for the uniform 3D limit Q122]) . 
These results show that the critical temperature T c (f2) never vanishes. In fact, as the rotation 
rate increases, the transition temperature becomes an increasing function of rotation rate, 
see Fig. [18] This is consistent with studies of BCS theory in the solid state setting beyond 
the semi-classical Gorkov approximation 225J. The increase in transition temperature at 
high rotation rate is related to the enhancement of the density of states at the Fermi level, 
due to the 1/y/E singularities of the quasi-lD motion along the rotation axis. These results 
show that the upper critical rotation rate found in Ref. 223 is a function of temperature; as 
T — » the BCS mean-field theory predicts that the system remains superfluid for all kpa s 
and Q. 

That BCS mean-field theory predicts the groundstate to be superfluid for all fcpa s and 
is a statement that the BCS state has a lower (free) energy than the normal phase for any 
attractive interaction. Despite the orbital motion induced by the rotation, the Fermi gas will 
always undergo a superconducting (SC) instability at sufficiently low temperature. This does 
not rule out the possibility of destruction of the superfluid phase by transition to another 
state. In Ref. [224] it was shown that the BCS phase competes with a charge density wave 
(CDW) state. The CDW is a state in which the atomic density is spontaneously modulated 
along the rotation axis with a period such that the atomic filling factor in each layer is 
v a = 2, see Fig. [19] This competition was analysed in detail for the case of interacting 
fermions in the LLL, allowing fluctuations of both CDW and SC. The calculations show 
that the dominant instability is towards the CDW state. (The transition temperature into 
the CDW is plotted as a dashed line in Fig.[TH]) Thus, an upper critical rotation frequency 
does exist on the BCS side of the transition: it is set by the transition into this CDW[224 . 



56 




Figure 18: Solid lines: Transition temperature for formation of a BCS superfluid state for 
a rotating Fermi gas in the uniform 3D limit, as a function of rotation rate f2 for different 
interaction strengths /cfOs- For weak interactions, as the rotation rate increases at fixed 
chemical potential /i (right to left on the graph), the transition temperature falls abruptly 
at the semi-classical upper critical rotation rate. However, it remains non-zero and rises as 
fl/jj, increases. For strong interactions the transition temperature becomes a monotonically 
increasing function of il. Dashed lines: Transition temperature for formation of a CDW 
state. [Reprinted figure from: G. Moller and N.R. Cooper, Phys. Rev. Lett 99, 190409 
(2007). Copyright (2007) by the American Physical Society. ] 




Figure 19: Schematic diagram of the charge density wave state that sets the upper critical 
rotation rate for destruction of fermionic superfluidity. The wavevector of the density wave 
lies along the rotation axis. It's period is set by the particle density, and is such that the 
atomic filling factor of each layer is v a = n a /{2ir£ 2 ) = 2. It is a fully gapped state. 
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At intermediate f2 it was argued that CDW and SC can co-exist, leading to a novel form of 
supersolidity. 

Transition from BEC to BCS 

The above considerations are expected to be valid on the BCS side of the transition. 
Understanding how the results connect to the BEC side is an interesting open problem. 
One thing that is clear is that, in the rapid rotation limit, in which all fermions (or bosonic 
molecules) are confined to the LLL, the transition from BCS to BEC side must involve a 
phase transitionF"! 

In the quasi-2D case this can be seen by noting that states on the two sides of the 
transition are topologically distinct. On the BCS side, where the atoms experience weak 
attractive interactions, the groundstate is a compact droplet with v a = 2 (a full Landau level 
for both atomic species); this state has two gapless edge modes, corresponding to density and 
spin excitations. On the BEC side, one anticipates that the groundstate is a Laughlin state of 
bosonic molecules, with v m = l /2| 21 l this state has only one gapless edge mode. The change 
in edge mode structure shows that the states on the two sides are topologically distinct; they 
must be separated by a (quantum) phase transition [226, 224, 227], An effective theory has 
been derived for this phase transition 227 , but the nature of the transition remains unclear. 

The existence of a phase transition in the 3D case was predicted [224\ based on the above 
2D argument combined with the observation that the 3D system forms a CDW of layers of 
v a = 2 states at high rotation rate. If the CDW order survives with the same period on 
passing to the BEC side, each layer must undergo the above topological phase transition; the 
only alternative is that the CDW does not survive (or changes period), which also signals a 
phase transition of another kind. 

6 Experimental Aspects 

Throughout the preceding text, existing experimental studies relevant to the theoretical 
results presented have been mentioned. In this section, I focus only on those aspects of 
experiment that relate to the strongly correlated phases at small filling factors. There also 
remain many interesting issues to explore experimentally beyond this strong-correlation 
regime. 

6.1 Experimental Signatures of Strongly Correlated Phases 

Experimental detection of the properties of cold atomic gases is typically performed following 
the removal of the trapping potentials and expansion of the gas. An important aspect of 
the physics of rapidly rotating atomic gases is that, provided the interactions are weak 
enough that the gas is in the LLL limit (|2~l|) , the expanded gas provides a direct map of the 
wavefunction prior to expansion: the only effect of expansion is a rotation and rescaling of 
the co-ordinates [228]. The use of expansion as a "wavefunction microscope" makes probing 
the rotating gas particularly straightforward. Contrast the case of a non-rotating gas where 
(in the simplest cases) the expanded image is a measure of the velocity distribution in the 
trap, and spatial structure is recovered only after Fourier transformation. 

The most dramatic feature of the strongly correlated phases predicted by theory is that 
they are incompressible. As a result, the density distribution in a trapped system will 
show characteristic "wedding cake" structures [125] . with regions in which the 2D particle 
density is pinned to values v*n v where v* are the filling factors at which the groundstate 
is incompressible. This density distribution differs markedly from the expected density 

20 This transition cannot be captured within the theory in the narrow resonance limit 12211 . as the bosonic 
molecules are non-interacting so the Laughlin state does not appear. 

21 The relation v m = u a /4 follows from the fact that there are half as many molecules as atoms, and the 
vortex density for a molecule is twice that for an atom [119] . 
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distribution of the vortex lattice phase which is compressible and hence follows a Thomas- 
Fermi distribution 43, 41]. An extreme case is the Laughlin state, v — 1/2, for which the 
average particle density is constant out to the radius R ~ 2\fNi. 

An interesting feature of particles in the 2D LLL is that the condensate fraction of 
the many-body phases can be determined from measurements of the expectation value of 
the particle density [229J . In this way, and by repeated measurements to construct the 
expectation value of the density, one could directly distinguish between a condensed state 
(such as the vortex lattice) and any strongly correlated phase for which the condensate 
fraction vanishes. 

In order to detect the nature of the correlated phase, it was proposed that the particle 
distributions in single-shot measurements of the expanded wavefunction be used to directly 
measure the density-density correlation functions of the states in the trap [228]. This method 
- also referred to as measurements of "noise correlations" 12301 - is particularly powerful in 
the present case, where the scaling of the wavefunction allows direct measurements of spatial 
correlations in the gas. The statistics of the density fluctuations can be used to measure 
the compressibility of the fluid, and to distinguish between different strongly correlated 
phases [225] , Noise correlations are also expected to show signatures of strongly correlated 
phases of bosons on rotating lattices [1541 [31] 

Further information on the nature of the incompressible phases can be deduced by study- 
ing the properties of their excitations. 

The clearest indication of the formation of an incompressible liquid phase in semicon- 
ductor systems is the appearance of a (fractionally) quantized Hall resistance|23j. This is 
a consequence of the pinning of the quasi-particles by disorder and the dissipationless flow 
of the incompressible liquid past these localized states. Proposals have been made [1581 [3T] 
for how analogous transport measurements might be performed in rotating Bose gases on 
optical lattices. 

Quantum Hall states, while gapped in the bulk, carry gapless edge excitations. These 
are collective modes (waves of particle density) on the surface of the incompressible fluid. 
The form of these excitations is fixed by the nature of the bulk phase [105] , Direct detection 
of the spectrum of edge excitations [23 lj could allow identification of the bulk incompressible 
state[i"U6]. 

A direct measurement of fractional statistics of the quasiparticle excitations of the in- 
compressible states would be a very dramatic confirmation of theoretical expectations. In- 
deed, even in conventional electronic FQHE systems, which are very well understood and 
much explored experimentally, an (accepted) experimental measurement of the fractional 
exchange statistics is still lacking. An ambitious scheme for measuring the semion statistics 
of the quasiparticle excitations of the Laughlin state of rotating bosons has been proposed 
in Ref.[235]. 

6.2 Strongly Correlated Phases: Experimental Status &; Outlook 

The reader is referred to Ref. [I] for a recent authoritative review of experiments on rotating 
atomic gases. Here I describe the current experimental status in connection with the strongly 
correlated phases which are predicted to occur at low filling factors, and discuss the outlook 
for future experimental work in this area. 

To date, the highest vortex densities in rotating atomic Bose gases have been achieved 
by the group of Eric Cornell. In these experiments, the mean interaction energy fj, ~ gn 
is sufficiently small that the gases are in or close to the 2D LLL regime, /i ;$ 2Tiuj±, Tvjj\\ . 
However, the smallest filling factor reported is v ~ 500 90 . This is deep in the regime v > v c 
where one expects the groundstate to be a vortex lattice [22l [49J, [50], and indeed this is what 
is observed. Reducing the filling factor further (reducing TV at fixed N v ) leads to a gas that 
is so dilute that collisions are too infrequent for thermalization on the lifetime of the gas. 
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Clearly further special efforts are require in order to reach the regime v < 6 where one 
expects quantum fluctuations of vortices to destabilize the lattice. There appear to be two 
main difficulties in reaching this regime: 

Weak Interactions 

In the 2D LLL, the typical interaction energy that stabilizes the strongly correlated 
phases is small. The relevant energy scale is the pseudo-potential Vo, Eqn. (f3"Tj) . To gain 
some physical intuition for the origin of this energy scale, it is useful to note that it sets 
the chemical potential (the mean interaction energy) for a rotating gas at a filling factor 
v ~ 1. The centrifugal limit (|19|) sets a maximum rotation rate = luj_, and hence a 
maximum vortex density (|10j) n v = 2Mto±/h. For v = 1 the maximum 2D particle density 
is 7i2d = n v = 2Mujj_/h and so the maximum 3D density is of order n ~ 2Mu>±/(han). The 
chemical potential associated with this density is 



h a s 1 a s 

V~ 9 n U — 2 huj ± ' m 

M a^a\\ an 



which recovers Vq pip up to numerical factors. 

This energy scale sets the size of the incompressibility gap of the strongly correlated 
phases. It therefore sets the required temperature scales in order to observe these phases. 
In typical experiments in magnetic traps, huj± ~ 5nK, while a s /a» ~ 0.01, leading to gaps 
that are small fractions of a nano-Kelvin, much smaller than the temperatures that can 
currently be reached in such experiment. 

Eqn. (|123[) highlights ways in which to overcome this difficulty: 

• Increase a s . This can be achieved by using a Feshbach resonance to enter a strongly 
interacting regime 233 . For bosons three-body losses can be large close to a Feshbach 
resonance, which may limit the usefulness of this approach. Losses are much smaller 
for Fermi gases close to an s-wave Feshbach resonance |234j . The FQHE regime for 
rotating fermions in the BCS to BEC crossover is therefore a very interesting situation 
to consider. It may be easier to achieve the strong-correlation regime (v ~ 1) for a 
strongly interacting two-component Fermi gas (close to the unitary limit) than for a 
typical Bose gas |235j . 

• Decrease ay . The introduction of an optical lattice, with wavevector parallel to the ro- 
tation axis can slice the condensate up into many parallel layers with a much stronger 
subband confinement (much smaller a[| ayer ), Fig. [T4l[126j . This can lead to a significant 
increase in the interaction energy scale |125] . It is important that the subband confine- 
ment retains rotational symmetry, to prevent loss of angular momentum or heating 
(see below). 

• Increase lo±. A tighter harmonic confinement would allow an increased vortex density 
and hence stronger interactions at v ~ 1. Much tighter traps can be formed by using 
atom chip traps formed on surface microstructures[236j, allowing transverse trapping 
frequencies on the order of 100kHz. (Such systems are highly anisotropic, wu <C lo_\_, so 
could potentially be used to study the 3D LLL regime |35|) .) Very tight confinement 
can also be achieved in the wells of an optical lattice. This might allow studies of 
rotating clusters of small numbers of atoms. There are interesting proposals 237|. [238] 
for how clusters of small numbers of particles might be transferred adiabatically from 
a non-rotating state into the strongly correlated Laughlin state. 



Optical lattice systems. The energy scales that set the correlation energies in optical 
lattices are large. Physical rotation of the optical lattice[134, 26 involves the same 
centrifugal limit (fTT)|) , again limiting the vortex density to the values leading to Q123P . 
This method will therefore likely also need to be supplemented with tight harmonic 
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confinement in order to achieve a large interaction energy at v ~ 1. However, the 
proposed schemes for imprinting Peierls phases onto the inter-site tunnelling matrix 
elements |137[ 128 ] [29 ] appear to allow field configurations that mimic uniform rotation 
without a centrifugal limit. If flux densities of order one flux per site can be generated, 
these systems could show strongly correlated states with a large energy scale (see S}3]) . 

Heating 

An atomic gas that has come to equilibrium at a rotation rate £1 is at rest in the rotating 
frame of reference. For example, in the vortex lattice phase the positions of the vortices 
are at rest in the rotating frame. In the laboratory frame, the fluid a distance rj_ from the 
rotation axis is moving at velocity of f2r_i_. As a result, the rotating gas will experience any 
non-axisymmetric perturbation that is static in the laboratory frame (e.g. a trap that is 
not circularly symmetric, or any scattered light from optical traps) as as a time-dependent 
potential in the rotating frame. This can lead to excitation, or heating, of the rotating gas. 

It is instructive to estimate the importance of such excitations on an incompressible 
liquid state, by applying a form of Landau criterion. For an incompressible liquid state, the 
collective excitations (in the bulk) have a gap, A, that is of order the interaction energy, Vo, 
or of order the cyclotron energy, 2hoj±, whichever is the smaller^] Therefore, the maximum 
possible gap (achievable for Vo ~ hui±) is A ~ Hlu±. The edge of a rotating cloud of radius 
R moves with a linear velocity of v = flR, which is close to v ~ cu±R for a system with a 
large number of vortices (|93|) . Due to this motion, in this edge region the collective modes 
with wavevectors q Z A/(hv) ~ 1/R have negative energy in the laboratory frame. These 
excitations can therefore be spontaneously generated by any perturbation that is static in 
the laboratory frame and has Fourier components on the relevant scale, q Z 1/R. This sets 
a strict requirement on the potential: any static perturbation in the laboratory frame that 
has a wavelength less than of order the radius of the cloud can efficiently generate collective 
modes the rotating gas (i.e. cause heating). 

The use of geometric phases to simulate rotation (or a uniform magnetic field) - either 
in a lattice-based system {T37| [28] [29] or in the continuum |239j - has the advantage that the 
gas remains at rest in the laboratory frame. This removes the relative motion of the gas 
with respect to possible perturbations at rest in the laboratory frame, and eliminates these 
as sources of heating. Experimental work has recently introduced a vector potential |240j of 
vanishing curl using two-photon induced geometric phases, setting the stage for extensions 
to vector potentials with non-zero curl. 

7 Summary 

Rapidly rotating degenerate atomic gases offer the possibility to explore the physics of 
quantized vortices and vortex lattices in novel parameter regimes. This review has described 
some of the novel equilibrium properties that can appear: for rotating spinless and multi- 
component Bose gases; for Bose gases on rotating lattices; and for rotating Fermi gases. 
Interesting new physics can emerge both within the realm where mean-field theories are 
applicable (as is dramatically illustrated by the novel vortex lattice phases in Figs. BJ [T5l 
and I16|). and within regimes where mean-field theory fails (as evidenced by the theoretical 
predictions of some very exotic strongly correlated phases). While theoretical studies have 
helped to clarify some of the expected properties of these systems, our understanding is 
still very much limited. It is hoped that these theoretical predictions will motivate and 
inform future experimental investigations of these novel phenomena. At the very least, 
these theoretical studies show that very interesting aspects of many-body quantum physics 

22 The lowest energy gap is proportional to Vo in the 2D LLL, Vo <C TujJ±. For Vo 2> hiOj_, there will be 
mixing of Landau levels and the lowest energy gap will not grow beyond the cyclotron gap 2hu>j_ . 
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are at play in rotating degenerate atomic gases; this is an area where future theory and 
experiment are likely to uncover dramatic new phenomena. 
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A Haldane pseudo-potentials 

Within the space of 2D LLL states, any two-body interaction that is rotationally symmetric 
in the 2D plane can be conveniently parameterized by its "Haldane pseudo-potentials" [211 . 
Owing to the rotational symmetry, the interaction conserves the relative angular momentum 
of two particles. Thus, its effects are determined by the expectation value of the two-body 
potential in the two-particle wavefunction 

* mc ,m(ri,r2) <x (C 1+ C2) mc (Ci-C 2 ) m e - |Cl|2/4 e - |C2|2/4 e ^ /2a " e - zi/2a » (124) 



Vrr, 



J /V(n-r 2 ) \y mc , m (ri,r 2 )\ 2 d 3 ri d 3 r 2 (125) 



where m and m c are the relative and centre-of-mass angular momenta. For spinless bosonic 
(fermionic) particles, m must be even (odd) so only the subset of even (odd) Haldane pseudo- 
potentials of the interaction are important. For contact interactions (|13[) . the only non-zero 
pseudo-potential is 

47r?i 2 a s f ,, 1 _qiv2 , 2w„2 1 t^, 2 12 h 2 a R 



M J (ttoI) 2 (7ra|) V Trajan v ' 

This is the only energy scale entering the low-energy physics of a rapidly rotating Bose gas 
in the 2D LLL regime. 

For the 3D LLL regime (|35|). an approximation often made [38] is that the particle density 
in the z direction is set by the Thomas- Fermi distribution. (See £12.61 for a discussion of this 
approximation.) One can view this as a change in the subband wavefunction in the z 
direction, ip(z). In the Thomas-Fermi approximation, the normalized particle density in the 
z direction is taken to be 



^ tfWI -mv-w m ) Qzl - Wz) (127) 

in which the extension along the z axis, W z , is set by the chemical potential, \i = (1/2)Mlo^W 2 . 
The only influence of the subband wavefunction in the z-direction ip(z) is to introduce a nu- 
merical factor to the Haldane pseudo-potentials, which for contact interactions is / dz\ip(z)\ 4 . 
Using the Thomas-Fermi form for the particle density in the z direction leads to 
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B Numerical Techniques 



Determining the groundstate of a rapidly rotating atomic gas in a harmonic well poses 
an essentially strongly interacting, non-perturbative theoretical problem. In the absence 
of interactions, the many-particle states at high angular momentum are degenerate, £)2.H 
Interactions lift this degeneracy and select a groundstate. To understand the properties of 
the groundstate (and the low-lying excitations) it has proved very useful to make use of 
exact diagonalization studies. 

The restriction of single particle states to the 2D LLL ([22]) limits the size of the Hilbert 
space. This makes exact diagonalization feasible for systems that can be large enough to 
observe clear signatures of strongly correlated many-particle states. Numerical studies have 
been performed in several different geometries, described in detail below. In each case, the 
strategy is the use of the "configuration interaction" method: 

• One chooses a finite number of (orthonormal) single particle basis states, <j> a (r), with 
a = 1, 2 . . .M which span the 2D LLL for the geometry considered. For example, in 
the uniform plane, these are the 2D LLL wavefunctions of Eqn. (l2"2l with a labelling 
the angular momentum rn = 0, 1, ... M — 1, where M sets the cut-off in total an- 
gular momentum. The wavefunctions in other geometries are discussed below. The 
interaction Hamiltonian may then be written in a second-quantized form as |241j : 

1 M 

V = - ^ V a p jS al e aj s a 1 as (129) 

a. /3, 7, 5—1 

where fiw are the field operators for particles (which could be bosons or fermions) in 
the single particle state (f> a (r). The matrix elements are 

V aPl& = (a0\V\Sj) = Jd 3 r J dV^(r)^(r')<MO<M»-) V(r - r') (130) 

where V(r) is the two-body interaction potential. 

• One generates a complete set of (orthogonal) many-particle basis states for a system of 
N particles occupying the single particle basis Q (r), with a = 1,2... M. These states 
are conveniently taken to be the Fock states, represented by the occupation numbers 
of the M orbitals, \rii, 112, ■ ■ -tim)- (The Fock states are symmetric/antisymmetric 
under particle exchange for the case of spinless bosons/fermions.) For bosons, the 
total number of such many-particle states is 

(N + M-l)\ N,M— too -yjy mi s 

AH(M-l)! ' { } 

where the limit N — > 00 is taken with fixed v — N/M, which sets the value 7. The 
number of basis states can be reduced by the use of symmetries, which depend on 
the geometry (see below). Nevertheless, the general feature that the number of basis 
states grows exponentially with system size survives, and sets the ultimate limitation 
to exact diagonalization studies. 

• Within the (symmetry reduced) many-particle basis, one constructs the matrix el- 
ements of the Hamiltonian (|129p . This is typically a very sparse matrix. Efficient 
numerical solution is possible for a small number of eigenvalues/eigenvectors using the 
Lanczos method. This allows the diagonalization of systems with of order 10 7 basis 
states on a modern desktop PC. 
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The different geometries that have been used in such studies are the "disc", "sphere" 
and "torus" geometries, and are described in detail in i]B.llB.2IB.3l They each have their 
own benefits and/or drawbacks. 

• Disc. This is appropriate for the direct simulation of experimental systems, in which 
particles are confined in a harmonic well. It allows studies of the effects of inhomo- 
geneity, including edge structures. Limitations arise in the study of correlated phases 
which can appear at large angular momentum. The inhomogeneity makes it difficult 
to observe emergent bulk phases. For a finite-size system, it is unclear how to define 
the number of vortices N v , so the filling factor (|67[) is not accurately defined. 

• Sphere. This is useful for studying bulk incompressible liquid phases. The uniformity 
of the geometry reduces the finite size effects associated with numerical studies of the 
inhomogeneous system on a disc. The number of particles N and vortices N v are fixed 
integers, allowing a well-defined filling factor v to be studied, via 



The offset S is referred to as the "shift" , which is a characteristic of a given incompress- 
ible liquid phase [232]. Owing to this offset (|132[) . different competing incompressible 
liquid states at a given filling factor can be studied independently on the sphere if they 
have different shifts. This is an advantage if the shift of the groundstate is known: 
consistent extrapolation to the thermodynamic limit along a sequence with that shift 
is strong evidence for this phase describing the thermodynamic limit. However, it is a 
disadvantage if the possible groundstates (and their shifts) are unknown. The sphere 
is not convenient for the study of crystalline states. 

• Torus. This is useful for studying bulk incompressible liquid phases, and bulk crys- 
talline phases (states with broken translational invariance). Finite size effects are 
reduced by the periodic boundary conditions. The number of particles N and vortices 
N v are integers, allowing a well-defined filling factor v to be studied 



Since there is no shift, all possible groundstates in the thermodynamic limit can com- 
pete within the same finite-size calculation. This geometry is therefore useful if the 
nature of the phase in the thermodynamic limit is unknown. For crystalline phases, 
the possibility to vary the geometry (aspect ratio) of the torus is a powerful way to 
investigate the optimal translational symmetry of the crystal. Certain incompressible 
liquid states are characterized by a groundstate degeneracy on the torus [ITS] : the ap- 
pearance of this non-trivial degeneracy in the exact spectrum is strong evidence for 
the groundstate being in this topological phase. It is somewhat harder to implement 
calculations on the torus than on the disk or sphere, especially when full advantage is 
made of the translational symmetries 243 . 

In the following I provide specific details for the different geometries that are required 
for the implementation of exact diagonalization studies. 



N v = -N-S. 



(132) 



v 



B.l Disk 



The orthononormal single particle basis states in the 2D LLL are 




(134) 
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with m = 0, 1, . . .. 

For contact interactions (|13[) . the matrix elements of the Hamiltonian are 

(mi + 777,4)! 

2mi+T7i2 



V - V v 1 I 4; A CI ^ 

"mi,m2,m3,m4 — *0 r,m,4-mo u m 1 +m 2 ,m 3 +m4 j \ ±00 J 



where Vq is the Haldane pseudo-potential (|31|) . 

For a set of iV particles, in orbitals rrii, the total angular momentum 



N 

l = y^TOj 

i=l 

is conserved by the interactions. The many-particle Hilbert space can therefore be split into 
subspaces of fixed L and N. There is an additional symmetry reflecting the independence of 
the interaction energy on the centre-of-mass co-ordinate of the particle. This symmetry can 
be used to further reduce the size of the Hilbert space, but this is not commonly implemented. 

B.2 Sphere 

A monopole of N v = 25 (N v integer) flux quantum is chosen to be positioned at the centre 
of a sphere, and to have uniform flux density over the surface [2 11] . The radius of the sphere 
is set by the condition that the total flux is 25 = AirR 2 n v , i.e. R = v~S£ where I is the 
conventional magnetic length (f25|) . The lowest-energy single particle states, analogous to 
the lowest Landau level states, consist of 25 + 1 degenerate levels. Using a gauge choice 
that is symmetric about the polar axis, the single particle wavefunctions on the sphere may 
be written [2H] 



25 + 1 / 25 



4-7T V S + rn 



1/2 

u s+m v s -" 1 (136) 



for m = -5, -5 + 1 . . . S, where u = cos(6»/2)e^/ 2 , v = sin(6y2)e^/ 2 and (0, (f>) are the 
polar angles on the sphere. 

The matrix elements of the contact interaction (|13j) are 

„ _ v [(25 + 1)!] 2 (25 + mi + m 2 )\{2S - mi - m 2 )\ 

'mi,m2,m3,m4 — ^0 i ^7711+7712, 7713+7714 ; I J 

5(45 + lyJutiiS + rm)\(S - mi)\ 



where we have chosen a 2D contact interaction of strength g 2 & = g/(v2iTai\) to take account 
of the quasi-2D motion along the rotation axis. 

These interactions (and any rotationally invariant interaction on the sphere) preserve 
the z-component of the angular momentum 



L,=^mi. (138) 

i=l 

This can be easily imposed to reduce the size of the many-body basis. Rotational invariancc 
also conserves the total angular momentum. The Hilbert space can therefore be further 
reduced to the many particle basis states of definite L 2 . In practice it is difficult to construct 
a basis of fixed L 2 , so this is not commonly used. 

B.3 Torus 

The "torus" geometry is defined by a unit cell on the plane on which opposite faces are iden- 
tified to impose periodic boundary conditions. It is specified by two (linearly independent) 
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basis vectors a, b in the plane. We shall focus on the simplest case of a rectangular unit 
cell, with basis vectors ax and by, with x, y orthonormal. A consistent quantum theory 
on the periodic geometry requires that an integer number of flux quanta N v pierce the unit 
cell, with N v = abn v — ab/(2ir£ 2 ) an integer. To construct single particle states it is con- 
venient to use the Landau gauge, with a vector potential directed along y. Then the linear 
momentum along y is conserved and spans the N v states in the 2D-LLL. Expressing the 
linear momentum along y in units of it can take the integer values m — 0,1, ..N v — 1. 
The single particle states in the 2D LLL arc 245J 



exp(— z 2 /2a|) 
(W) 1/4 



1 



1/2 



E 



exp 



. (X m + pa)y (X m + pa- xf 



£ 2 



2£ 2 



(139) 



where X m = 2n£ 2 m/b. 

The matrix elements of the contact interaction are 



V, 



mi ,?7i2 jWls .m^ 



v ^Vb E < 

P,q=— oo 



7 mi+m2 ,m3+m4 



(140) 



where S^J = 1 if i — j = modiV v and otherwise. 
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The total momentum in the y direction 



K v = E mi 



(141) 



is conserved modulo N v . This provides a convenient way in which to reduce the size of 
the many-body basis size (by a factor order N v ). Owing to the symmetry under magnetic 
translations along x, a second conserved momentum K x can be constructed [243]. At a filling 
factor N/N v — p/q (with p and q co-prime), with number of particles N = pN and number 
of flux iV v = qN, the momenta (K x ,K y ) take N 2 distinct states in a (rectangular) Brillouin 
zone. Construction of many-particle basis states that are eigenstates of both K x and K y 
allows a significant further reduction in the overall size of the Hilbcrt space. At particular 
points in the Brillouin zone there exist additional point-group symmetries. 
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